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I.  INTRODUCTION 


It  is  generally  known  that  a periodic  array  of  linear  scatterers 
(dipoles)  may  yield  a unity  reflection  coefficient  at  a certain 
frequency  (resonance),  while  the  complementary  configuration  consisting 
of  slots  yields  a unity  transmission  coefficient  at  the  same  fre- 
quency [1,2].  This  statement  is  based  on  Babinet's  principle,  i.e., 
it  is  true  only  when  the  two  complementary  arrays  are  infinitely  thin 
and  are  made  of  perfectly  conducting  material.  If,  however,  the 
structures  are  of  finite  thickness  or  consist  of  more  than  one  layer 
of  dipoles  or  slots,  Babinet's  principle  no  longer  applies  and  the 
reflected  signal  from  the  dipole  configuration  no  longer  equals  the 
transmitted  signal  of  the  slot  configuration.  The  effect  of  screen 
thickness  on  a single  slotted  surface  has  been  investigated  earlier  [3,4] 
where  it  was  shown  that  the  main  effect  of  increasing  the  screen  thick- 
ness is  to  make  the  resonance  curve  narrower.  In  another  report  the 
reflection  characteristics  from  an  arbitrary  number  of  layers  of' 
resonant  dipoles  (loaded  as  well  as  unloaded)  is  found  and  it  is 
demonstrated  how  band  filter  characteristics  can  be  obtained  by  proper 
design  [5].  Finally  the  transmission  properties  for  two  slot  arrays 
separated  by  a dielectric  slab  has  been  obtained  [6].  It  indicates 
that  a band  filter  curve  can  be  obtained  but  that  certain  problems 
occur  for  high  angles  of  incidence  (80°)  in  the  slot  H-plane.  It  was 
felt  that  an  improved  band  filter  curve  with  a flatter  top  could  be 
obtained  by  the  addition  of  dielectric  layers  on  the  outside  of  the 
biplanar  slot  configuration.  This  report  considers  this  problem, 
and  it  will  be  shown  that  the  expectations  were  indeed  true.  The 
effects  of  surface  waves  on  the  outside  dielectric  layers  are  that  of 
producing  nulls  in  the  transmission  curve,  and  it  will  be  shown  how  to 
locate  these  nulls  outside  the  passband.  Similarly  it  will  be  shown 
how  a related  phenomenon  in  the  middle  layer  of  dielectric  can  produce 
a null  in  the  transmission  curve.  This  is,  however,  not  a true  surface 
wave  but  can  be  avoided  in  the  same  manner. 


II.  SOLUTION 


A.  Calculation  of  the  Induced  Voltages 


Consider  Fig.  1 showing  two  slotted  arrays  mounted  behind  each 
other.  Array  No.  1 is  located  in  the  XZ-plane,  while  array  No.  2 is 
separated  from  array  No.  1 by  a dielectric  slab  of  thickness  d2  and 
relative  dielectric  constant  e2*  Further,  a dielectric  slab  of 
thickness  di  and  relative  dielectric  constant  ei  is  mounted  in  front 
of  array  No.  1 and  a slab  d3(e3)  is  mounted  behind  array  No.  2. 


r 
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DIELECTRIC 


Figure  1.  General  view  of  a bi -planar  slot  array  sandwiched 
between  three  dielectric  slabs  d-j , d2  and  d^ 


Each  array  contains  (2R+1)  rows  and  (2K+1)  columns  of  slots  each 
of  length  2t  and  with  interelement  spacings  Dx  and  Dz.  The  slots  in 
arrays  No.  1 and  No.  2 are  loaded  with  the  load  admittances  Y[_]  and  Yl2» 
respectively.  A plane  wave  is  incident  upon  this  configuration  at  an 
angle  <f>j  measured  from  the  negative  Y-axis  in  the  XY-plane  (E-plane 
or  4>-plane)  or  at  an  angle  e-j  measured  from  the  negative  Y-axis  in  the 
YZ-plane  (H-plane  or  e-plane).  We  seek  the  field  transmitted  through 
this  configuration. 

We  now  denote  the  vector  effective  height  of  the  reference  slot, 

No.  00,  in  array  No.  1 for  h|*(e.f).  This  effective  height  has  been 
determined  in  Append ix_A.  The  induced  current  in  this  reference 
element  is  then  h(e-j).Hi»  where  Hi  is  the  magnetic  field  vector  of  the 
incident  field.  Since  array  No.  2 is  shielded  from  the  incident 


2 


field  by  array  No.  1,  no  current  will  be  directly  induced  in  array 
No.  2 due  to  the  incident  field.  If  we  further  recall  the  definition 
of  mutual  admittance  between  slots  we  can  thus  readily  write  the 
following  two  equations  for  the  reference  slots  in  array  No.  1 and 
array  No.  2,  respectively: 


(1) 


■ <\l  + VA> 


t VT  v<2> 

i2  oo 


l 


(2)  0 - yJ,  »<;>  + (YL2  * yJJ)  V<2> 


where 


= mutual  admittance  sum  between  the  reference  element  oo 


nm 


in  array  n and  all  the  elements  in  array  m,  i.e., 


(3) 


T R K T 

<nm  = l l n rk  e' 

r=-R  k=-K  n,rK 


(4) 


r R K r 

yA=  l l Vorke 
M r=-R  k=-K  0,rK 


where 


(5) 


eDx  r s in<f>. 


$ = < 


k sine.. 


/ (n ) 


for  4>-pl ane  scan 
for  e-plane  scan 


and  further  Vv7  denotes  the  terminal  voltage  across  element  rk  in 
array  n,  wher^  because  of  Floquet's  theorem 

(6)  •<&>.-»  • 


The  mutual  admittance  sum  Y-j^  and  the  admittance  sum  of  a single  slot 
array  coated  with  dielectric  and  backed  by  a ground  plane  (the 
superscript  G)  has  been  derived  and  discussed  in  detail  in  another 
report  [8]. 

Equations  (1)  and  (2)  formally  determine  the  unknown  quantities 
VAq  and  V^'.  However,  since  we  are  at  present  only  interested  in 
the  transmitted  field  determined  entirely  by  the  voltages  in  the 
second  array,  we  shall  determine  only  those  voltages.  Fromtqs,  (1) 
and  (2) : 


■—  — 


_ 


(11)  Fe2  * KSiT  C’  - cos^fte  - Fels1n6Iel- 


The  transmission  coefficient  for  the  above  biplanar,  dielectric 
coated  slot  configuration  is  now  defined  as  the  ratio  between  HF-S- 
as  given  by  Eq.  (8)  above  and  the  field  Hfq.Qp.  transmitted  through  the 
"Equivalent  Opening"  defined  as  an  aperture’with  the  area 


Thus,  for  the  transmission  coefficient  T we  find  by  division  of  Eqs, 
(13)  by  (8) 


1/T  = 


r i K 

U1- 1 -jgd^cos  j0 
BZ  A H . cos)  . [ e * ^ 

0 1 I 


2N  FE3,H3[Fe)  -vf  i> 


Substituting  Eqs.  (7)  and  (12)  into  Eq.  (14): 


(15)  1/T  = 


6Zo  Dx  Dz  c0S[e.J  e 


,1  -36d2  cos|6J 


2 h?<e1>Ps<VFE3,H3|Fel  -rf  Fe2 


fVfll+VLl  ^VA3+VL2-'‘V12V21 


The  vector  effective  height  ht^e.)  of  a dielectric  covered  slot 
has  been  derived  in  Appendix  A as 


_n  „ COS  BA l-  COSB«.  n 

(A-16)  - - e,  I FE1H1  pjle,) 


where  p5(e ) is  the  radiation  pattern  of  the  individual  element  under 
transmitting  conditions. 


Substituting  Eq.  (A-16)  into  Eq.  (15)  yields 


JK,  K?  Z" 

(16)  1/1  ' I ~ 4"  Y^dl  ,2,3;YL1 ,YL2* 

rEl  ,H1  E3,H3 


5 


' ' - * - _ 


....  S 


where 

(17) 


■n[a/\)‘ 


C0SgAJ>-C0Ses,e1 

F Yl2  F 

L singze  j 

Lei  Fe2J 

yi  DJi  cos^l 

(18)  7K  2 = D,  . IT  7 1 

PS^ei)  Pt(et) 


(!9)  Y^dl  ,2,3’  YL1’YL2 


tYAl+YLl^YA2+YL23  - Y12Y21 


T VT  -jgcLcc 


As  shown  in  Appendix  B,  Eq.  (B-22) 


(B-22) 


4 F 


,n<-FP 


El, HI  FE3,H3  2 G^(0)1/2  G^(0)1/2  E1 ’H1  ~ E3’h 


where  (0 ) and  G^O)  are  defined  below. 

Substituting  Eq.  (B-22)  into  Eq.  (16)  yields  the  very  simple 
formula : 


(20)  i it  - Y(dl,2,3;YLl’YL2)  P - 

n ' 2 G^<0)V2  Gfi.(0)l/2  <-FEl,Hl<-FE3,H3  • 

The  admittances  in  the  expression  for  Y(d]  2 3iY|l»YL2)  given  by 
Eq.  (20)  above  has  been  derived  earlier  as  [8] 


(21) 


'A1 


= I 


gai (n) 


+ J B 


A1 


(22) 


A2 


= I G^(n) 


+ J B 


A2 


where 


pG 

GA1 


(n)  = 


the  conductance  of  a propagating  mode  of  slot  array 
No.  1 coated  with  a dielectric  layer  of  thickness  d] 
and  dielectric  constant  e]  and  backed  by  a ground 
plane.  In  particular,  n=0  corresponds  to  the 
principal  (desired)  propagation,  while  other  values 
of  n (if  any)  correspond  to  grating  lobes. 


B 


^ = the  total  susceptance  of  a slot  array  coated  with  a 

dielectric  layer  of  thickness  d]  and  dielectric  constant 
C1  and  backed  by  a ground  plane  at  a distance  d2  and 
where  the  "cavity"  is  filled  with  a relative  dielectric 
constant 


G G 

G/\2  and  B/\?  are  defined  in  an  anologous  way  for  array  No.  2.  Note, 


however,  that  array  No.  2 is  coated  with  the  dielectric  d^Ug). 

Finally  the  mutual  admittances  y]^  and  Y^  have  been  determined  in 
Appendix  C,  and  found  to  be  entirely  imaginary:  Yjo  = j Q-j 7 » where  Q-j 2 
is  given  by  Eq.  (C-7).  By  substituting  Eqs . (21),  (22)  ana  (C-6)  into 
Eq.  (20)  we  obtain 


(23) 


1/T  = 


QGA1  Cn>+j (bAi+VLi  ) ]CZGA2C n)+j (BA2+YL2^  ^+(^l  2^21 


_ n 


y a 

-jed2c°s1e 

l 1 J <-F.,  ,,, <-Fr 


El, HI  E3,H3 


Equation  (23)  may  be  greatly  simplified.  If  the  two  arrays  are 


identical,  but  Yi -j ^Y«  2 » it  can  be  shown  that  we  always  have  a lossy 
jfficient  [6 " 


transmission  coefficient  [6].  If  the  two  arrays  as  well  as  the  load 
admittances  Y[_i  and  Y|_2  are  different  it  is  not  definite  that  loss 
will  always  result.  However,  there  seem  to  be  no  special  advantages 
in  this  case  and  it  will  not  be  pursued  further  at  this  time. 


Yl_rYL2' 


We  shall  finally  consider  the  Synmetric  Case:  d-j  =d3 , e]  03. 


I 


In  the  symmetric  case  Eq.  (23)  above  reduces  to 


1/T  = 


+ J(BA+YL)32  + Q2  -j6d2cos(  0 J 


2j  Q Gjj(0) 


rEl ,H1  * 


In  order  to  find  the  extrema  of  Eq.  (24)  we  now  find  the  numerical 
value 


(25)  VlTr  = ~~c  — 2 o 
4G^(0)2Q2  L 

where  for  brevity  we  have  put 


B4  + 2B2 


^lG®(n)]2  - QZ)  + ( [iG^(n)]2  * Q 2f 


6 - B“+Yl. 


Inspection  of  Eq.  (25)  shows  that  it  contains  three  variables: 

B,  zGWn)  and  the  coupling  Q.  Of  these,  the  parameter  B will  vary  by 
far  the  most  as  a function  cf  frequency  except  when  we  are  operating 
close  to  a grating  lobe  which  will  be  investigated  separately.  Thus, 
in  order  to  determine  the  extrema  of  1 / j T 1 2 given  by  Eq.  (25),  it  is 
permissible  to  differentiate  with  respect  to  B (provided  Q and  G?(n) 
vary  slightly: 

(27)  dtumiL , (KGG(n)]v) 

A 


which  has  the  roots 


B0/Q  - o 


[^(n)]2 


8 


- 1 Jit1  .*■  .LflgSgiig11-* 


The  transmission  coefficients  T0  corresponding  to  the  root  b0,  and 
T+1  corresponding  to  the  roots  B+-| , are  obtained  from  Eq.  (25)  as 

1 Iga(")2  n 1 

(30)  l/|TJ=i  ~4 + — p5— 

° 2 GG(0)Q  GG(0)J 


iGfl(n) 

(31)  1/|T+1I=-| 

- GG(0)  . 

Thus,  from  Eq.  (31)  above  we  observe  that  a unit  trans- 
mission coefficient  can  be  obtained  only  at  T+^  for 


(32)  GG(0)=£GG(n) 


i.e.,  if  no  free  space  grating  lobes  exist.  Apparently  this  condition 
is  independent  of  the  coupling  Q.  However,  note  that  B+,  exists  only 
for  Q>EG)((n).  (This  makes  physical  sense  since  grating'lobes  means 
losing  energy  which  must  be  taken  from  the  principal  direction.) 

Similarly  we  observe  from  Eq.  (30)  that  the  transmission  coef- 
ficient T0  will  in  general  have  values  smaller  than  unity,  except  for 

(33)  Q=lGG(n), 

which,  by  inspection  of  Eq.  (29),  is  seen  to  correspond  to  B+,=0,  i.e., 
the  three  roots  B+^  and  Bq  will  in  that  case  merge  into  just~Ane  root. 

Finally,  from  Eq.  (29)  we  observe  that  if  Q <G^(n),  B+,  will  simply 
not  exist,  while  B will  yield  a lossy  transmission  coefficient  cf. 

Eq.  (30).  0 

A pictorial  summary  of  the  findings  above  is  shown  in  Fig.  2. 

This  clearly  shows  that  we  will  obtain: 


0 > G« (0) 
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III.  MORE  DETAILS  CONCERNING  THE  TRANSMISSION  COEFFICIENT 


The  investigation  above  has  been  conducted  in  terms  of  the 
parameters  G^(n),  B and  the  coupling  Q.  The  first  two  were  investi- 
gated in  great  detail  earlier,  while  the  coupling  Q is  investigated 
in  Appendix  C.  Based  on  these  findings  it  is  clear  that  the  specific 
values  B+,  and  B0  can  be  obtained  for  several  frequencies  leading  to 
rather  complicated  transmission  coefficient  curves.  More  insight  into 
this  complex  problem  is  obtained  by  studying  Fig.  3 (typical  for 
({.-plane  scan)  and  later  Fig.  4 (e-plane  scan).  Here  curve  "a"  depicts 
a typical  performance  of  B as  a function  of  frequency.  It  is  observed 
that  B has,  at  least,  two  poles: 

1.  At  the  onset  of  the  first  grating  lobe  in  the  dielectric  e2  between 
the  two  slot  arrays. 

2.  At  the  frequency  where  a surface  wave  propagates  along  the  slab 
d-j  (or  dg)  with  dielectric  constants  (or  eg). 

Since  e2,  for  reasons  to  become  clear  later,  usually  is  larger  than 
ei  , the  onset  of  the  first  grating  lobe  in  e2  will  occur  at  a lower 
frequency  than  the  onset  of  the  first  surface  wave  in  the  outer 
dielectric  layers  d]  and  d3-  It  is  also  observed  that  B attains  the 
value  B+,  and  B several  times  as  a function  of  frequency  as  stated 
above. 

Further,  curve  "b"  in  Fig.  3 depicts  the  coupling  Q as  a function 
of  frequency.  As  shown  in  Appendix  C,  the  leading  term  in  this  coupling, 
valid  for  no  grating  lobes  in  dg,  is  given  by 

A V \ 


Pre2(02^ 


/IJsin^Jl  ^DxDz 


SinUgdgSggJ 


A close  examination  of  this  term  will  show  that  it  varies 
relatively  slowly  as  a function  of  frequency,  in  particular  for 
02d2s32'VlT/2 . It  is  also  seen  to  be  negative  imaginary.  However,  when 
approaching  a grating  lobe  condition  in  the  middle  slab  d2,  a positive 
pole  occurring  at  the  onset  of  the  grating  lobe  in  the  ({.-plane  will 
completely  dominate  the  leading  term  referred  to  above  (Eq.  (34)). 

Also  at  some  frequency,  before  the  actual  onset  of  grating  lobe,  it 
will  cancel  that  term  making  the  coupling  Q equal  to  zero. 

Finally  Fig.  3c  shows  a typical  transmission  curve.  It  is 
observed  that  a unity  transmission  coefficient  is  obtained  every  time 
B attains  the  values  Bf  and  a somewhat  lower  value  is  obtained  at  B . 
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Note  also  that  zero  transmission  is  obtained  when  the  coupling  Q attains 
the  value  zero  (Luebbers  anomaly)  and  when  a surface  wave  is  excited 
in  the  outer  dielectric  slabs.  However,  at  the  onset  of  grating  lobes 
in  the  middle  slab  d2»  we  observe  a singularity  in  B as  well  as  the 
coupling  Q.  The  fact  that  both  of  these  occur  at  the  same  frequency 
can  be  shown  not  to  cause  any  discontinuity  or  singularity  in  the 
transmission  coefficient  as  seen  in  Fig.  3c. 

A closer  investigation  will  show  that  the  first  (lowest  frequency) 
hump  of  the  double  hump  resonance  can  be  designed  to  be  the  most  stable 
as  a function  of  the  angle  of  incidence. 

Similarly  Fig.  4 shows  trie  susceptance  B and  the  coupling  Q 
typical  for  e-plane  scan.  Note,  that  no  discontinuity  occurs  at  the 
onset  of  grating  lobes  on  the  middle  layer  (d2)  (because  of  the  pattern 
factor  P2(02)'  anci  alsc  that  surface  waves  will  occur  at  a much  higher 
frequency  than  for  <j>-plane  scan  shown  in  Fig.  3 (thus,  they  are  not  shown 
in  Fig.  4).  However,  at  the  onset  of  a free  space  grating  lobe  we 
observe  a pole  in  B but  not  in  Q.  Consequently  a null  will  be  observed 
in  the  transmission  curve  |Tj  as  shown.  This  is  quite  remarkable  since 
the  onset  of  a grating  lobe  in  the  e-plane,  without  an  outer  dielectric 
layer  d"| , produces  no  null  but  only  a discontinuity  in  the  derivative 
of  | T | . However,  in  the  4>-pl ane  a null  (Woods  anomaly)  will  be 
observed  without  an  outer  dielectric  layer,  but  as  seen  in  Fig.  3 above, 
with  a dielectric  we  observe  only  a discontinuity  in  the  derivative  of 
|T[.  In  other  words,  the  dielectric  can  displace  the  classical  Woods 
anomaly  from  the  ^-piane  to  the  9 -pi ane.  Beyond  this,  the  outside 
dielectric  layer  is  a media  for  possible  surface  waves  in  particularly 
in  the  6-plane  as  was  shown  earlier. 
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Figure  4.  Typical  behavior  of  the  parameters  pertinent  to  the  trans- 
mission coefficient  for  the  e-plane  (H-plane): 
a.  The  total  susceptance  B;  b.  The  coupling  Q between 
the  two  slot  arrays;  c.  The  resulting  transmission 
coefficient  |T|. 


IV.  HOW  TO  DESIGN  A TWO  LAYER  BAND  PASS  FILTER 


From  the  previous  discussion  we  should  now  be  able  to  design  a two 
layer  band  pass  filter  with  a transmission  curve  shaped  as  shown  in 
Fig.  2 (and  more  specifically  in  Figs.  3 and  4).  This  design  will  be 
characterized  essentially  by  two  frequencies  with  unit  transmission 
coefficient  and  one  frequency  in  between  where  the  transmission  coeffi- 
cient is  down  slightly  below  unity.  Beyond  these  resonances  the  surface 
waves  in  the  dielectric  will  give  rise  to  a number  of  resonances  as 
discussed  above,  however,  they  will  usually  be  so  narrow  and  change  so 
fast  with  angle  of  incidence  that  their  significance  in  practice  is 
quite  small . 

A.  The  Inter  Element  Spacings  Dx  and  Dz 


It  is  generally  true  that  the  larger  Dx  and  Dz,  the  narrower  the 
band  pass  curve  and  vice  versa  for  smaller  interelement  spacings. 
However,  the  onset  of  the  grating  lobes  is  also  associated  with  Dx  and 
Dz.  Since  this  phenomenon  usually  leads  to  nulls  or  spurious  resonance 
close  to  the  onset  frequency  (in  particular  when  dielectric  layers  are 
involved),  we  usually  have  certain  constraints  placed  on  Dx  and  Dz  by 
the  grating  lobe  considerations . 

In  order  to  obtain  the  proper  coupling  between  the  two  arrays, 
the  dielectric  constant  e2  between  the  two  slot  arrays  will  generally 
be  chosen  higher  than  the  dielectric  constants  e]  and  e3  in  the  outer 
layers,  as  explained  later.  Therefore,  the  onset  of  grating  lobes  will 
occur  first  between  the  arrays  as  we  go  up  in  frequency  rather  than  in 
the  outer  layers.  Furthermore,  as  explained  above  and  also  in  Appendix 
C,  the  onset  of  grating  lobes  between  the  two  arrays  is  preceeded  by  a 
null  in  the  coupling  between  the  two  slot  arrays  resulting  in  a null 
in  the  transmission  curves  (Luebber’s  anomaly).  This  anomaly  occurs 
at  a frequency  approximately  5-8%  below  the  onset  of  the  first  grating 
lobes  in  d2,  when  the  onset  frequency  is  determined  by  the  well  known 
grating  lobe  formula 


(35) 


3-10 


gr 


Dx(sinPi  + JT) 


Thus,  if  we  decide  to  place  the  Luebbers  anomaly  higher  than  the 
upper  band-pass  frequency  f , we  must  require 

u p 


(36) 


D <*v 

x 


3-10 


1.1  f (sin  ,,  *J c2) 
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The  absolute  lower  limit  for  Dx  (and  Dz)  is  of  course  observed 
when  the  elements  interfere  physically  with  each  other.  However,  at 
least  two  factors  might  prevent  us  from  pushing  to  this  limit: 

a)  The  bandwidth  of  the  transmission  curve  might  become  too 
broad, 

b)  The  transmission  coefficient  at  higher  frequencies  might 
become  undesirably  high  due  to  the  fact  that  higher  order 
resonances  in  the  elements  are  not  attenuated  properly  by 
the  presence  of  grating  lobes,  and 

c)  mechanical  weakness  of  the  screen. 

Thus,  in  conclusion,  we  can  state  that  Dx,  in  most  practical  cases, 
should  be  chosen  approximately  10-15%  smaller  than  the  right  hand  side 
of  inequality  (Eq.  (36)).  If  transmission  phase  equality  between  the 
two  principal  planes  is  important,  Dx  and  Dz  should  be  chosen  such 
that  the  amplitude  of  the  transmission  curves  for  the  two  planes  are 
as  much  alike  as  possible  as  explained  in  Section  5.  The  condition  of 
(Eq.  (36))  still  holds  but  the  best  values  Dx  and  Dz  will  often,  in  this 
case,  have  to  be  finally  determined  by  calculation  of  the  complex 
transmission  curve  as  given  by  Eq.  (24). 

B.  Choice  of  c-j  (eg)  and  d^  (dg)  for  the  Outer  Layer 

One  of  the  fundamental  problems  of  space  filters  made  from  layers 
of  periodic  surfaces  is  the  fact  that  the  bandwidth  changes  dramatically 
with  angle  of  incidence.  This  can  be  traced  to  the  fact  that  the  real 
part  of  the  scan  admittance  changes  with  angle  of  incidence.  As  pointed 
out  in  an  earlier  report  [8],  one  way  to  make  this  scan  admittance  constant 
with  scan  angle  is  to  place  a dielectric  slab  on  the  outside  of  the  slot 
array.  When  placed  directly  on  top  of  the  slots,  the  thickness  of  the 
slab  that  will  yield  the  most  constant  admittance  was  determined  to  be 
somewhat  larger  than  A]/4  and  the  value  of  the  dielectric  constant  should 
be  chosen  ei=1.2-1.8.  Values  of  d]  and  e]  in  that  range  will  provide 
good  compensation  for  the  scan  admittance  for  scan  up  to  as  high  as  1 80° 
in  both  the  <f>-  and  0-plane. 

C.  Choice  of  £2  and  dg  for  the  Middle  Layer 

While  the  outer  dielectric  layer  plays  the  role  of  producing  a 
constant  bandwidth,  the  role  of  the  middle  layer  is  to  produce  the  proper 
coupling  between  the  two  slot  arrays.  As  explained  above,  if  the 
coupling  is  too  small,  we  will  observe  undercritical  coupling  with  loss. 

If  the  coupling  is  too  strong,  overcritical  coupling  with  two  resonance 
frequencies  with  too  deep  a null  in  between  will  result.  In  both  cases, 
the  loss  is  given  by  Eq.  (30)  which  has  been  plotted  in  Fig.  5 as  a 
unction  of  Q/Gx(0).  By  inspecting  Eq.  (30)  we  further  note,  that 
values  of  Q/Gx(0)  > 1 correspond  to  the  loss  observed  between  the  two 


maxima  for  overcritical  coupling  while  values  of  Q/G(?(0)  < 1 correspond 
to  the  loss  at  the  single  resonance  for  undercritical  coupling.  In 
most  applications  we  would  strive  for  a design  ranging  from  critical 
to  slightly  overcritical  coupling  with_ typically  0.5  dB  loss  cor- 
responding to  Q/G^(0)=1  and  Q/Gj((0  )=%/2 , respectively  (see  Fig.  5). 

Ip  order  to  facilitate  the  design  we  have  in  Fig.  6a  through  6f  shown 
G^(0)/G^(0)frge  Spare  as  a function  of  di/q  for  various  angles  of 
incidence  on  bottv  planes  for  e-j  =1 . 05  through  2.0.  Inspection  of 
Eq.  (20)  in  Reference  [8]  shows  that  G^(0)/G^(0)free  Space  has 
extrema  forr . 't 
1+1 1 

Woskf=  *12- 


This  explains  why  the  extrema  in  Figs.  6 occurs  at  d] / At  =0 . 25  for  normal 
angle  of  incidence  but  at  higher  values  of  d^ / Aq  for  oblique  angles  of 
incidence.  Similarly  we  show  in  Fig.  7a  through  7j  the  coupling 
Q/G^(0)free  spaCp  as  a function  of  d£/ ^2  various  angles  of  incidence 
and  for  several  fixed  values  at  e 2-  Further,  it  should  be  noted  that  whil 
G^(0)/G^(0)free  spaCe  is  independent  of  the  interelement  spacings  Dx 
and  Dz,  this  is  nbt  the  case  for  Q/Gg(0)fr  space'  In  fact,  certain 
values  of  Dx/A  and  Dz/a  can  at  certain  angles  of  incidence  produce  zero 
coupling  as  seen  for  example  in  Fig.  7h  for  <J> -}  =80°  resulting  in  zero 
transmission,  see  Fig.  3,  (Luebber's  anomaly).  Consequently  the  design 
curves  for  Q/G^(0)f  space  have  been  Pl°tted  for  various  typical 

values  of  Dx/A  = Dz/a.  From  Eq.  (C7)  we  observed  that  the  leading  term 
for  Q/Gj>(0)free  space  is  proportional  to 


when  Dx/A  = Dz/A  is  small,  i.e.,  the  extrema  for  normal  angle  of 
incidence  is  obtained  for  d2/A2=0.25,  while  it  moves  to  higher  values 
of  d2/>2  f°r  oblique  angles  of  incidence.  However,  for  larger  values 
of  Dx/A  = Dz/A  some  increased  spread  in  Q/G^(0)free  space  W1th  angle 
of  incidence  is  observed.  This  spread  will  prove  quite  useful  in 
obtaining  an  optimum  filter  design. 

n 

The  use  of  the  design  curves  for  G(0)/G^(0)free  space  and 
Q/Gn(0)free  cpaqe  now  9°es  ^ll<e  this:  Let  us  assume  that  we  wish 

to  design  a band  filter  with  critical  to  overcritical  coupling  where 
the  dip  in  the  middle  is  down  no  more  than  0.5  dB.  From  Fig.  5 we  see 
that  the  ratio  Q/G^(0)  should  be  between  one  (critical  coupling)  and 
1.4  (overcritical  coupling).  In  order  to  obtain  constant  bandwidth  we 
must  first  of  all  choose  e-j  between  1.2  and  1.8  as  discussed  above. 

We  next  must  choose  a coupling  Q/G^(0)free  sp  from  the  design 
curves  Fig.  7a  through  7j,  which  has  a variation  with  angle  of  inci- 
dence matching  that  of  G^(0)/G^(0)free  s such  that  their  value  is 
between  unity  and  1.4.  ■ 


resonance  for  undercritical  coupling. 


figure  7.  The  normalized  coupling  Q as  a function  of  the  electrical 
thickness  d2/*2  of  the  middle  dielectric  slab  for  various 
angles  of  incidence.  Dielectric  constants  e 2 and  inter- 
element spacings  Dx  = Dz  as  follows: 
a:  = 1 .00  D.  = D_  = 0.300a. 
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Figure  7.  The  normalized  coupling  Q as  a function  of  the  electrical 
thickness  d2/A2  of  the  middle  dielectric  slab  for  various 
angles  of  incidence.  Dielectric  constants  e2  and  inter- 
element  spacings  Dx  = Dz  as  follows: 
e:  e2  = 1 .90  Dx  = Dz  = 0.320A. 
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Figure  7.  The  normalized  coupling  Q as  a function  of  the  electrical 
thickness  d2/X2  of  the  middle  dielectric  slab  for  various 
angles  of  incidence.  Dielectric  constants  e2  and  inter- 
element spacings  Dx  = Dz  as  follows: 
h:  to  * 1.9  D = D_  = 0.394A. 
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Figure  7.  The  normalized  coupling  Q as  a function  of  the  electrical 
thickness  °f  the  middle  dielectric  slab  for  various 

angles  of  incidence.  Dielectric  constants  c2  and  inters 
element  s pacings  Dx  = Dz  as  follows: 
i:  e2  = 2.3  Dx  = D7  = 0.320a. 
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Figure  7.  The  normalized  coupling  Q as  a function  of  the  electrical 
thickness  d2/A2  of  the  middle  dielectric  slab  for  various 
angles  of  incidence.  Dielectric  constants  e2  and  inter- 
element spacings  Dx  = Dz  as  follows: 
j:  e0  = 2.3  D = D = 0.355A. 


V.  PHASE  DELAY 


So  far  we  have  concentrated  our  attention  almost  entirely  on  the 
amplitude  of  the  transmitted  signal  both  as  a function  of  frequency 
as  well  as  angle  of  incidence.  However,  for  many  applications,  as  for 
example  a radome,  the  variation  of  the  phase  with  angle  of  incidence 
is  of  importance.  From  Eq.  (24)  the  complex  transmission  coefficient 
T+i  at  the  two  frequencies  f+1  when  unit  transmission  is  always 
obtained:  " r r ' u 


g ;i°) . 


ga<°> 


• 71  , 

J 2 - d^cos 
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As  discussed  in  detail  above,  the  condition  for  critical 
coupling  yielding  a flat  top  is  given  simply  by 

(38)  |Q|  = Gj((o). 

This  condition  (Eq.  (38))  makes  the  three  frequencies  f+,  and  f 
degenerate  into  each  other  yielding  0 

I h-ll 


1/T+  = e 

-1  crit 
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Inspection  of  Eq.  (37)  shows  that  the  phase  delay  at  the  resonant 
frequency  can  be  broken  into  several  components.  First,  we  have  the 
exponential  factor  where  the  7r/2-term  can  be  traced  back  to  the  coupling 
between  the  two  slot  arrays  while  the 

f^cos^g  -term 

k i. 

j 

is  caused  simply  by  the  physical  relocation  of  the  second  ("transmitting") 
array  with  respect  to  the  first  ("receiving")  array.  Note  that  this 
factor  is  apparently  independent  of  polarization  (assuming  critical 
coupling  for  the  two  polarizations) . Second,  in  Eq.  (37)  we  have 
twice  the  phase  of  the  factoi  £1  ,H1 • This  factor  was  simply  defined 
as  the  complex  factor  by  which  the  incident  field  should  be  multiplied 
in  order  to  account  for  the  dielectric  slab  c] , di . Similarly  the  field 
transmitted  was  multiplied  by  F^3,h3  due  to  the  slab  £3,  d3-  In  our 
situation  the  space  filter  is  made  symmetrically,  i.e.,  e-|  = £3,  d]  = d3 
which  explains  the  factor  of  two.  It  is  now  pertinent  to  observe  the 
phase  of  F£i  m as  a function  of  angle  of  incidence  for  various  slab 
thicknesses  d^  as  well  as  dielectric  constants  as  can  be  seen  in 
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Fig.  8.  We  observe  that  <-F^-|  1n  general  will  differ  in  the  two 
polarization  planes.  However, ’if  for  example  e]  = 1.30  and  d-| / A-j  =0,4, 
we  obtain  very  nearly  the  same  phase  variation  in  the  two  planes  for 
the  angle  of  incidence  up  to  around  70°. 

One  may  be  tempted  to  conclude  from  the  above  analyses  that  use 
of  such  a thickness  is  all  that  is  required  to  insure  the  same  phase 
delay  for  the  two  polarizations.  This  is,  however,  not  the  case.  A 
third  factor  must  be  considered,  namely  the  variation  of  the  resonant 
frequency  f^-,  =f 0 with  angle  of  incidence.  This  point  is  perhaps  best 
understood  if  we  consider  the  phase  variation  of  a typical  critical 
coupled  space  filter  as  shown  in  Fig.  9.  Within  the  passband  the 
phase  is  seen  to  vary  almost  linearly  as  a function  of  frequency. 

However,  at  the  resonant  frequency  f±i=f0  the  phase  delay  has  been 
given  by  Eq.  (37)  leading  to  the  same  value  for  the  two  polarizations 
provided  <F^i  = <F(-|l . Thus,  as  seen  from  Fig.  9,  a difference  in  the 

resonant  frequencies  f^  and  foe  leads  to  a difference  in  phase 
between  the  two  planes.  One  might  agree  that  a difference  in  resonant 
frequency  for  the  two  polarizations  could  be  compensated  for  by  a similar 
but  opposite  difference  between  <Ff\  and  <Fhi . While  there  seems  to 
be  no  physical  reason  to  prevent  the  application  of  such  a scheme,  we 
have  so  far  not  been  able  to  take  full  advantage  of  such  an  approach. 

We  assume  this  is  because  the  variation  of  the  resonant  frequency  with 
angle  of  incidence  is  very  complex  and  depends  on  the  outer  as  well  as 
the  middle  dielectric  layer  and  the  interelement  spacings  Dx  and  Dz. 

We  may  conclude  that  the  design  procedure  outlined  above,  con- 
cerning e],  d-|  and  e2,  d£,  is  necessary  to  obtain  a good  band  filter 
curve  at  all  angles  of  incidence,  but  it  is  net  quite  sufficient  to 
insure  a constant  resonant  frequency  with  the  angle  of  incidence,  which 
in  turn  is  necessary  to  produce  the  same  phase  delay  for  the  two 
polarization  cases.  We  shall  in  the  next  section  present  results  which 
show  how  additional  manipulation  of  Dx  and  Dz  can  yield  a good  compromise 
between  phase  and  amplitude  requirements. 
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Figure  8 (cont.) 
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figure  9. 


Typical  amplitude  and  phase  response  of  a bi -planar 
slot  configuration  illustrating  the  difference  in 
phase  between  the  4>-  and  e-plane  because  of  their 
difference  in  resonance  frequency. 
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VI.  RESULTS 


Based  on  the  design  criteria  outlined  above,  we  have  built  a 
dielectric  clad  bi-planar  slot  configuration  as  shown  in  Fig.  1. 

A computer  program  based  on  Eq.  (23)  has  been  written  that  yields 
the  magnitude  and  phase  of  the  transmission  coefficient  and  plots  the 
magnitude  in  dB  as  a function  of  frequency  for  various  angles  of  in- 
cidence in  the  E-  as  well  as  the  H-plane.  A copy  of  this  program  is 
given  in  Appendix  D.  A set  of  transmission  curves  of  a typical  design, 
denoted  P-27,  in  the  frequency  range  5 -18  GHz  for  various  angle  of 
incidence  is  shown  in  Fig.  10.  Figure  11  similarly  shows  an  expanded 
view  (8-12  GHz)  of  amplitude  as  well  as  transmission  phase  of  the  same 
design.  Measured  curves  obtained  by  our  swept  frequency  set  described 
earlier  [10]  are  shown  in  Fig.  12.  For  dimensions  for  the  P-27  design, 
see  Table  I and  Fig.  17.  It  is  observed  that  the  bandwidth  remains 
almost  constant  with  angle  of  incidence  all  the  way  up  to  80°  for  both 
the  principal  planes.  Recall  that  without  dielectric  compensation  the 
bandwidth  would  vary  approximately  like  l/cos280°  - 1:33(1).  We  also 
observe  the  nulls  in  the  transmission  curve  for  the  E-plane  U-plane) 
first  at  the  lower  frequencies  because  of  coupling  break  down  between 
the  two  arrays  and  at  higher  frequencies  because  of  the  surface  wave 
in  the  outer  dielectric  layers  as  discussed  in  Section  3.  In  between 
these  nulls,  resonances  with  unity  transmission  will  occur  as  discussed 
earlier,  however,  as  can  be  seen  from  Fig.  10,  they  are  very  narrow  and 
change  dramatically  with  angle  of  incidence.  In  practice  this  makes 
them  look  smeared  and  low  in  value  as  seen  in  the  measurements  in  Fig. 

12.  It  is  further  observed  that  the  transmission  curve  for  P-27  shown 
in  Fig.  11  has  an  undesirable  deep  valley  between  the  two  peaks  for  80° 
angle  of  incidence  in  the  E-plane  (<j>-plane).  As  discussed  earlier  (see 
Section  4.3),  this  is  caused  by  too  strong  coupling  between  the  two 
arrays.  By  increasing  the  interelement  spacings  Dx  = Dz  from  0.96  cm 
for  the  P-27  design  to  1.13  cm,  a new  design  called  P-24  is  obtained. 

(For  dimensions,  see  Table  I and  Fig.  17.)  By  comparing  Fig.  7e  and 
Fig.  7g  depicting  the  coupling  Q between  the  two  arrays  for  Dx  = Dz 
= 0.96  cm  and  1.12  cm,  respectively,  it  is  observed  that  such  an  increase 
in  interelement  spacings  will  result  in  a lower  coupling  Q,  in  particular 
at  high  angles  of  incidence.  The  calculated  transmission  curves  for  this 
new  design  P-24  is  shown  in  Fig.  13  for  the  frequency  range  5- 18  GHz. 
Figure  14  shows  an  expanded  view  of  both  the  amplitude  and  phase  response 
for  the  band  pass  range.  The  measured  transmission  curves  of  P-24  are 
shown  in  Fig.  15.  It  is  observed  in  the  calculated  transmission  curves 
(Figs.  13  and  14)  that  the  deep  valley  at  80°  angle  of  incidence  in  the 
E-plane  is  completely  gone,  as  expected,  however,  the  measured  curve 
for  80°  angle  of  incidence  show  "wavey"  response  in  the  H-plane.  We 
attribute  this  to  measuring  problems  peculiar  to  high  angle  of  incidence. 
It  is  also  observed  that  the  first  nulls  (Luebbers  anomaly)  in  the 
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E-plane  curves  have  moved  down  considerably  in  frequency  as  discussed 
in  design  constraint  for  Dx  and  Dz  in  Section  4.1.  If  the  amplitude 
response  is  the  only  consideration  for  a potential  application  of  any 
of  these  two  designs,  the  P-24  design  in  Fig.  14  may  well  be  deemed 
superior  to  the  P-27  design  in  Fig.  11.  However,  if  the  panel  is  going 
to  be  used  for  building  a radome,  for  example,  the  phases  of  the  trans- 
mitted signal  must  also  be  composed  within  the  band-pass  region  of  the 
two  designs.  Thus,  the  phase  response  of  P-27  is  shown  in  Fig.  11  and 
of  P-24  in  Fig.  14.  It  is  seen  that  for  a fixed  frequency,  the  P-27 
design  shows  less  phase  variation  with  angle  of  incidence  and  polarization 
than  does  P-24.  For  the  sake  of  comparison,  we  show  in  Fig.  16 
the  amplitude  and  phase  response  for  a typical  x/2-radome  (er  =4.2, 
no  loss).  Comparison  between  Fig.  11  for  the  P-27  design  and  Fig.  16 
for  the  x/2-radome  design  shows  that  although  the  phase  response  for 
the  two  polarizations  has  a crossover  point  close  to  the  resonance 
frequency  in  the  X/2-  case,  the  two  designs  are  otherwise  quite 
comparable.  If  the  interelement  spacings  Dx  = Dz  are  further  increased 
to  1.355  cm,  a third  design  denoted  P-8  is  obtained  (for  dimensions 
see  Table  I and  Fig.  17).  The  calculated  curves  are  shown  in  Fig.  18 
and  the  measured  in  Fig.  19.  It  is  clearly  observed  that  the  first  null 
(Luebbers  anomaly)  has  now  moved  down  into  the  pass  band  making  this 
design  undesirable  at  the  higher  angles  of  incidence  in  the  E-plane 
U-plane).  However,  such  a configuration  could  conceivably  be  used  as 
a filter  letting  the  low  angle  of  incidence  through  (0-30°),  but  stopping 
signals  arising  in  the  E-plane  at  angles  of  incidence  in  the  range  60-80°. 


Design 

Parameters 

TABLE  I 

for  Various  Filter  Designs 

P-27 

P-24 

P-8 

Dx  = Dz  (cm) 

0.96 

1 .13' 

1.355 

D]  = Di  (cm) 

1 .10 

1.11 

0 85 

D2  (cm) 

0.60 

0.72 

0.70 

cl  = e3 

1.30 

1.30 

1.50 

e2 

1.90 

1 .90 

1.90 

i (cm) 

0.34 

0,34 

0.375 

w ( cm) 

0.18 

0,18 

0.18 

T (cm) 

0.0071 

0.0071 

0,0071 

LTL  (cm) 

0.32 

0.32 

0.32 

ZTL  (n) 

240 

240 

240 

BDF 

1.26 

1 .26 

1,26 

b (cm) 

0.060 

0.060 

0.051 

c (cm) 

0.310 

0.310 

0.320 

d (cm) 

0.170 

0.170 

0.183 

Comparison  between  the  calculated  and  measured  curves  above  shows 
in  general  good  agreement  with  respect  to  resonance  frequency,  surface 
wave  nulls  and  coupling  nulls  (Luebbers  anomaly).  As  far  as  bandwidth 
is  concerned,  however,  the  calculated  curves  are  consistently  more 
narrow  banded  than  the  measured.  While  this  has  been  observed  at 
other  occasions  and  been  attributed  to  tolerance  problems,  we  believe 
another  essential  reason  is  the  fact  that  the  load  admittance  4/Z£  Z[_ 

has  been  treated  as  pure  susceptance.  While  this  has  been  a good 
approximation  in  the  past,  it  will  be  seen  by  inspection  of  a single 
slot  drawn  to  scale  in  Fig.  17,  that  the  ends  of  this  elements  center 
sections  are  quite  comparable  in  size  to  the  actual  element.  Since  the 
fields  at  the  ends  of  the  loads  are  even  stronger  than  in  the  actual 
slots,  it  is  obvious  that  by  no  means  insignificant  radiation  will 
take  place  at  the  end  of  the  loads.  This  has  a broadening  effect 
on  the  transmission  curve.  Although  this  discrepancy  could  be  accounted 
for  at  present,  it  will  be  much  simpler  to  let  the  new  improved 
Poisson  approach  automatically  incorporate  any  such  effect. 

In  the  Ku-band  region,  it  is  also  observed  that  the  high  calculated 
spurious  resonant  frequencies  on  the  c^-plane  in  the  measured  version 
has  the  character  of  a smear.  This  is  precisely  what  we  desire  in  many 
applications.  Further,  in  the  <t>-plane,  we  observe  some  measured  resonances 
in  Ku-band  not  seen  in  the  calculated  version.  This  is  simply  due  to  the 
fact  that  an  odd  mode  is  excited  which  presently  is  not  included  in  the 
solution.  Again,  it  will  be  seen  in  the  improved  Poisson  version. 
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Figure  13 


Design  P-24.  Calculated  transmission  curves  as  a function 
of  frequency  for  various  angles  of  incidence. 

Top:  E-plane  ( -plane).  Bottom:  H-plane  (e-plane). 
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Figure  15.  Design  P-24.  Measured  transmission  curves  as  a 
function  of  frequency  for  various  angles  of 
incidence. 

Top:  E-plane  U-plane). 

Bottom:  H-plane  (e-plane). 
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function  of  frequency  for  various  angles  of 
incidence. 

Top:  E-plane  (<f>-plane). 

Bottom:  H-plane  (e-plane). 
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Figure  16.  Amplitude  and  phase  response  for  a typical  A/2  radome 
Top:  Parallel  polarization  (E-plane  = 4>-plane). 

Bottom:  Perpendicular  polarization  (H-plane  e-plane) 
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VII.  CONCLUSIONS 


In  this  report  we  have  investigated  the  transmission  properties 
of  a metallic  radome  configuration  consisting  of  two  resonant  slotted 
conducting  surfaces  which  are  surrounded  and  separated  by  layers  of 
dielectric.  Only  the  symmetrical  case  consisting  of  two  identical  slot 
arrays  was  considered. 

In  an  earlier  report  it  was  found  that  a bi -planar  arrangement  of 
identical  slot  arrays,  without  dielectric  layers,  would  always  have 
some  transmission  loss  even  at  the  resonant  frequency.  However,  in 
this  report  it  was  shown  that  the  inclusion  of  appropriate  dielectric 
layers  to  form  a 5-layer  symmetrical  configuration  would  yield  a 
surface  with  good  bandfilter  characteristics.  That  is,  it  would  have 
a transmission  coefficient  curve  with  no  loss,  a relatively  flat  top, 
and  sharp  skirts  at  the  resonant  frequency  for  angles  of  incidence  up 
to  and  including  80°  in  both  the  E-plane  and  the  H-plane. 

The  phase  of  the  energy  transmitted  through  such  an  arrangement 
was  also  investigated.  It  was  found  that  by  proper  design  the  phase 
difference  between  TE  and  TNI  polarization  could  be  as  low  as  15°  for 
angles  of  incidence  up  to  60°  off-normal  and  over  most  of  the  passband. 

It  was  further  discovered  that  the  dielectric  layers  were 
instrumental  in  achieving  constant  bandwidth  as  a function  of  incidence 
angle.  The  thickness  of  the  layers  and  the  relative  dielectric 
constants  of  the  materials  was  important.  Typical  values  were  1.2  to 
1.8  for  the  relative  dielectric  constant;  total  sandwich  thicknesses 
generally  turned  out  to  be  slightly  more  than  a quarter  wavelength. 


APPENDIX  A 


VECTOR  EFFECTIVE  HEIGHT  OF  A 
DIELECTRIC  COVERED  SLOT 


_q  In  this  appendix  we  shall  determine  the  vector  effective  height 
hs(0i)  of  a slot  covered  with  an  infinitely  large  dielectric  slab  of 
thickness  d and  relative  dielectric  constant  c.  As  is  well  known 
such  a slot  is  the  equivalent  of  two  magnetic  dipoles,  one  on  each 
side  of  the  conducting  screen,  with  the  magnetic  currents 


(Al)  V(z)  = VsineuUe-|z|) 

as  shown  in  Fig.  Al . 

We  have  here  assumed  that  the  slot  is  transmitting  and  reasonably 
thin  such  that  a sinusoidal  voltage  distribution  is  a good  approximation. 

Further  is  the  propagation  constant  along  the  slot  as 
determined  earlier  [11,8].  In  the  same  references  we^have  also 
determined  that  the  far  field  from  a voltage  element  z dz  is  given  by 


where 


v<2)  h^rr  C0S8i  fh,e 


ie/2  -ft  coset  cos(sd  cosei  " 0/2) 
”J  e cose . sinU^d  coset) 

JT  cos*.  cos(ed  cos*.  - */2) 

eJ<^/  ^ - 1 I-,- L 

J cos*  sin(e  cl  cos*t) 


for  H-plane 


for  E-plane 


where 


6£  8 
sine- 


and  the  angle  e is  defined  by  the  equation 


/e  cose. 

(A6)  tan(gd  cose.  - e/2)  = — — ^ cot(g  d cose.) 

• cosy • 0 t 

and  4>  by 


(A7) 


tan(gd  cos^.  - 4>/2) 


fz  cos<j>. 

COS<j>t 


cot(g£d 


cos4>t ) . 


Thus,  the  far  field  for  a magnetic  dipole  located  immediately  in  front 
of  an  electric  screen  is  then  obtained  by  substituting  Eq.  (Al) 
into  Eq.  (A2)  and  integrating  from  to  +£: 


(A8) 


H = e.  ju)£Q  2 V ^ cose.. 


rSi  je  z sine.  n 

e e t sineD(^-|z|  )dz 
-i  e 


where  the  factor  of  two  in  Eq.  (A8)  is  due  to  the  contribution  from 
the  image  in  the  electric  screen. 


Tedious,  but  straightforward,  evaluation  of  Eq.  (A8)  yields  by 
noting  that  g^sine^.  = gsine.  : 


(A9) 


_ » i 4V  Fr  u COS0. 

H = e.  J-  1 

1 n„  .. 

. 2 

sin  8. 


/ D \2 

(f-)- 


-jer0 

e 

4nr 

o 


j j-[cos  g°AX.  cos  ( gas i ne . ) -cos  gD5.g  ] 


- sine.  singuA«.  sinles-sine^ ) 


The  transverse  component  Nt  of  N is  now  defined  by  [12] 


J ti  4Trrrt 

(A,0)  "t  * -r  Tk  * 

e 


Substituting  Eq.  (A9)  into  Eq.  ( AT 0 ) 


. 4V  F cose.  D n 

(All)  N.  = - e.  — — *-= r-  f-  [coseUAtcos(etsin0. ) 


r r5in  ei 


- cosB°te]-sinei  sineDAt  sin(etsinei ) . 


The  vector  effective  height  h“(e,)  of  a dielectric  covered  slot 


is  now  defined  by 


-0  Nt  Nt 

(Ai2)  ■ tt  - rpr 

in  V sine  t 


Substituting  Eq.  (All)  into  Eq.  (M2): 


n * i 2 F F iiCose. 

(»13)  H°s(6f)  - - 0,  i 


I 


j-  [coseDAt  cos(gtsinei.)  - co$eDte]  - sine .sing^tsinUtsine^) 


D D 

A , cose  At -cose  n 

■-6^2Fe.h-7 *-p? 

*-•»  A 


P“(ef) 


where 


(A14)  pJCe.)  = 


L) 

[|-  [cosgDA£cos(e£sine.)-coseD£e]-sineisineDA2sin(8£sinei )] 
is  the  normalized  pattern  function  excluding  Fc  u. 

t ,rt 

It  is  now  simple  to  show  that  for  gPl 


(A15) 


D D 

COS3  At-COSe  l COS BAt-COS et 

£.  /-v  e 


r sirA 


si  net 


if  the  correct  t is  used! 
e 


By  application  of  Eq.  (A15)  we  may  write  Eq,  (A13)  as: 


_n  „ i coseAt-coset  n 

(A16)  h>,)  - - e,  i 2F£iH  su|  P?(6f). 


(Bl) 

(B2) 


(B3) 


(B4) 


where 


— — _ 


APPENDIX  B 

Proof  of  [FE>H|2  = Rg/RA 


We  have  earlier  found  [11]  (for  H-plane  scan) 
nO  -je  2dcose. 

= Re  1^  Up  e £ t 

RA  1+p  -jg  2dcose 

1-p  e e 1 


Ize. 

1+p 


cos0t 

cose^ 


cose..  - Jc  cos 0^. 

p = 

cose^  + JT  coso^ 

Substituting  Eqs.  (B2)  and  (B3)  into  Eq.  (Bl)  yields 

Rn  cose.  cos0.+/Tcos0  +(cos0.-/7cos0.)e 

j£=te^5SiT  ~ — 5 — 

cosQj+Jecose^-icosQ.-Jzcose. )e 


-jg  2dcose 

e 

-je  2dcos0 

e 


t 


t 


_C°s0t 


= Re  ^ 


. cose.  cos(e£dcos0t)+jje^iiTsin(eedcoset) 

i - ^ ' 


cose^. 


cose 


U i 

ie  c-os-6-  cos(6cdcos0t)+jsin(eedcos0t) 


2 / \ 2 
cose  Jcos  ^*™\Hc5SS7  )*'n  ^edc°set) 

' coseT  )^=f= ~-^=- ~ cos(B-a) 

/ /COS0i  \ 2 2 

x/e  Jcos  (eedcos0t)+sin  (Mcoset) 
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where 


1 1 COS0 . 

C0S3  = g- cosCe^dcose^)  and  sing  = g-yr  — S0’~  sin(Bedcosot) 


1 COS0  , 1 

cosa  = Av/1"  coseT  cos(3edcoset)  and  stna  = ^ sin(g£dcoset) 


p f COS0 , V p 1 

B = [cos  (gedcoset)+e|— ^ j sin  (eedcoset)] 


2 

COS  0 j.  p p -j 

(B8)  A = [e  — 5 — - cosc(g  dcosej+sin  (3  dcosej] 
cos^e.  e 1 e t 

Substituting  Eqs . (B5)  and  (B6)  into  the  formula 
cos(g-oc)  = cosgcosa+singsina 

yields 

1 cose 

(B9)  cosU-a)  - m p ^ . 


Substituting  Eq.  (B9)  into  Eq.  (B4)  and  making  use  of  Eqs.  (B7)  and 
(B8)  yields 


(BIO) 


R° . 1 rs6t  f 


ra  a2  Vcosei  J ' 


We  have  earlier  found  [13] 


(Bll) 


je/2  coset  cos($dcose.  - e/2) 

Fu  3 ^ A f .*  n / A rt  \ 


cose.  sin(g  dcose J 

i e t 
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where 


cose. 

(B12)  tan(gdcose.  - e/2)  = -IF  ■■■■--  cot(g  dcosej  . 

i y cose.  e x. 


Substituting  Eq.  (B12)  into  the  formula 


cos  p = 


l+tan2p 


yields 


(B13)  cos($dcose.  - e/2)  = 


/cose.  \ ? 

1+<\wsii7  )cot  (sedcoset> 


Substituting  Eq.  (B13)  into  Eq.  (Bll)  yields 


fB14l  F - i ej0/2  1 

(B  4)  fh  J e cosei 


Comparison  between  Eqs.  (B14)  and  (BIO)  now  readily  yields 


(B15) 


F |2  =-A 
ra 


Equation  (B15)  can  be  shown  also  to  hold  for  E-plane  scan  such  that  we 
have  in  general 


(B16) 


R° 

if  |2.!A 

E,H  ra 
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F 


I 

[ ! 


Equation  (B16)  above  could  also  be  proven  much  simpler  by  simple 
energy  considerations.  Let  the  induced  current  in  the  reference 
element  for  an  array  without  dielectric  be  denoted  by  I and  similarly 
be  denoted  Iu  for  the  same  array  with  a dielectric  layer.  Since  the 
array  has  the  same  aperture  with  and  without  dielectric,  the  gain 
must  be  the  same  in  the  two  cases,  i.e.,  the  energy  received  in  the 
reference  element  must  be  the  same: 


r ■ 


( B 1 7 ) 


= ml 

ga 


However,  we  also  have  defined 


(B18) 


F = i— 
rE  ,H  I 


Substituting  Eq.  (B17)  into  Eq.  (B18)  yields 


(B19) 


,2  GA(n=0) 

FE,hI  = GA(n=0) 


Note  that 


■ i 


(B20)  GA(n=0)  = \ ^ Ra 

Zo 


because  GA(n=0)  radiated  only  to  one  side) 


and  as  shown  earlier 


<B2,)  - 

Substituting  Eqs . (B20)  and  (B21)  into  Eq.  (B19)  yields 

z?  fiJC  , 

(B22)  / . 1 

IFE,hI  2GA(n=0)  . 
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APPENDIX  C 

THE  MUTUAL  COUPLING  yL  BETWEEN  TWO  PARALLEL  SLOT  ARRAYS 


In  tnis  appendix  we  shall  determine  the  mutual  admittance  y|-| 
between  the  reference  slot  in  array  No.  2 and  all  the  other  slots  in 
array  No.  1.  As  usual  we  shall  make  use  of  the  fact  that  a slot  can 
be  represented  as  two  magnetic  dipoles  mounted  on  each  side  of  a 
perfectly  conducting  ground  plane.  In  order  to  make  things  a bit  more 
clear  we  have  in  Fig.  Cl  located  the  dipoles  in  array  No.  1 a distance 
a-j  in  front  of  the  ground  plane  No.  1 and  similarly  the  reference 
dipole  in  array  No.  2 a distance  a2  in  front  of  ground  plane  No.  2, 
where  later  we  shall  let  a^  and  a2+0  to  simulate  a true  slot.  Only 
the  dipoles  facing  toward  each  other  will  be  considered  since  no 
coupling  can  exist  between  the  external  dipoles  because  the  presence 
of  the  ground  planes. 

To  find  the 
[1]  the  voltages 
reference  dipole 

the  elements,  four  infinite  series  of  images  are  created  as  shown  in 
Fig.  Cl.  The  distances  between  the  reference  dipole  [2]  and  all 
these  images  are  shown  in  Fig.  Cl,  as  well  as  their  relative  amplitudes, 
which  are  readily  evaluated  when  it  is  recalled  that  the  reflection 
coefficient  for  a magnetic  dipole  image  in  an  electrically  conducting 
ground  plane  is  +1.  By  inspection  of  Fig.  Cl  and  by  recalling  the 
definition  of  mutual  admittance  we  can  readily  write 


mutual  admittance  Y^-j  we  now  impose  upon  the  dipoles 
Vv  ' and  seek  the  current  ij^'  induced  in  the 
ofKarray  [2],  Apart  from  tne  direct  coupling  between 


(Cl) 
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oo 
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a1 , a2  -*0 


co 

l [Y21  ((2n+l)d2-a1+a2) 
n=0  +Y21  ((2n+l)d2+ai+a2) 
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+Y21  ((2n+l)d2-a]-a2)] 
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5d-a,  + a2- 
-5d  + a,  + a2 


Figure  Cl.  Relative  location  of  amplitude  of  the  images 
of  two  magnetic  dipoles  at  1 and  2. 


’52 


sin*2  + n5  r 


7b 


In  Eq.  (Cl)  we  now  let  a-j  and  a^  and  obtain 

(C3)  Y?l  = 4 l Y21((2n+l)d?). 

n=0  c 

Substituting  Eq.  (C2)  into  Eq.  (C3)  yields: 

T 8 Y o , 


8¥e2 

1 

. 2 

sin  e2£ 

4dxd2 

I I 

pre2*e; 

s 

‘1  "3 


+ jl  l Pre2(e2)  £ ^B^n+l  )d2s42 


n-j  n4  42  n=0 


_ j ^ £ Pim2^e2^  y e“^2^2n+1 )d2s52 

n2  n5  s52 

Applying  the  surrmation  formula  for  a infinite  geometrical  series  with 
ratio  smaller  than  unity  (assuming  for  the  moment  that  media  2 is 
slightly  lossy  to  avoid  mathematical  embarrassment  (!))  we  readily 
find  for  the  summation  over  n: 


J _ '^2_  1 

21  ~ 2 2 

siV-e2*  g|DxDz 


l l 


pre2^e2^ 


n1  n3  32 


■2  2 32  ^ 


pre2^e2^ 


+ V y rre2x"2'  1 

n,  n4  s42  ^^4? 


I l 


pim2(e2) 


n2  n5  52 


'2252 ' 


_ ...... 


SfcM# 


or 


r 


Yi,  = J Q; 


where 


21  sin2^  g2DxDz 


-I  I 

- n]  n-: 


Pre2^e2^ 


sin(62d2s32) 


I I 


Pre2^02) 


sinh(62d2s42) 


pf m2  ^ e2  ^ 


s i nh ( 2 ^ 


Note  from  Eq.  (C5)  or  Eq.  (C7)  that  y|i  is  purely  imaginary  which 
agrees  with  physical  reasoning  since  no  energy  can  be  lost  between 
the  two  ground  planes.  Also  note  that  for  increasing  d2,  the  last 
two  double  summations  in  Eqs.  (C5)  or  (C7)  will  be  small  compared  to 
the  first  double  summation  because  sinh(g2d2s42)  and  sinh(82d2s52)_><,0> 
provided  that  S42  or  S52  does  not  assume  values  close  to  zero.  This 
last  condition  will  be  observed  at  the  onset  of  a grating  lobe  which 
will  make  infinite.  However,  because  the  first  (and  third)  double 
summation  is  negative  (fpr  0<82d2s32  ) while  the  second  is  always 

positive,  there  will  be  a frequency  before  onset  of  grating  lobe 
in  the  <t>-plane,  where  we  will  obtain  a null  of  Yl. . This  null  will 
be  seen  in  the  analysis  sections  to  create  a null  in  the  transmission 
coefficient.  It  was  first  calculated  by  R.O.  Luebbers  [10],  and  has 
since  been  called  a "Luebbers  anomaly". 

A physical  explanation  of  this  phenomenon  goes  like  this:  For 

spacings  tt/4<B2c*2s32<3tt/4  the  first  double  summation  in  Eq.  (C7) 
(consisting  of  only  one  term  for  no  grating  lobe)  is  dominating  pro- 
ducing an  inductive  coupling  between  the  two  slot  arrays  (-j  yields 
inductance  for  susceptances) . The  second  double  summation  represents 
the  stored  capacitive  energy  caused  by  the  "strip"  structure  of  the 
elements  while  the  third  double  summation  represents  the  stored 
inductive  energy  caused  by  the  "cutting  of  the  strips  into  slots".  In 
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genera]  both  of  these  summations  will  be  dominated  by  the  first 
summation  representing  the  stored  energy  between  the  two  arrays  but  at 
the  onset  of  grating  lobe  the  second  or  third  double  summation  becomes 
dominating  resulting  in  a null  as  explained  above. 


y 

l 
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APPENDIX  D 
COMPUTER  PROGRAM 


This  appendix  presents  Fortran  listings  of  the  main  program  and 
subroutines  which  have  been  developed  to  calculate  the  magnitude  and 
phase  of  the  transmission  coefficient  for  the  biplanar  slot  array 
sandwich  which  has  been  described  in  the  text.  The  main  program  is 
explained  in  detail  while  the  subroutines  are  only  briefly  described. 

A detailed  description  of  the  subroutines  appears  in  Reference  [8], 

A.  Main  Program 

A Fortran  listing  of  the  main  computer  program  which  calculates 
the  magnitude  and  phase  of  the  transmission  coefficient  is  given  at 
the  end  of  this  lection.  The  program  is  made  up  of  three  basic 
sections.  The  first  (lines  1-101)  sets  the  scan  plane,  the  frequency, 
and  the  angle  of  incidence  for  which  the  transmission  coefficient  will 
be  determined.  The  second  section  (lines  102-190)  calculates  the 
admittance  of  each  array  as  well  as  the  mutual  admittance  between  them. 
The  third  and  final  section  (lines  191-288)  finds  a total  admittance 
for  the  entire  sandwich  and  normalizes  it  to  yield  a transmission 
coefficient.  This  section  also  contains  provisions  for  printing  out 
the  magnitude  and  phase  of  the  transmission  coefficient  as  well  as 
plotting  the  magnitude  (lines  262-286). 

Now  considering  the  first  section,  lines  12-13  are  commands 
which  are  unique  to  the  system  for  which  the  program  was  written  and 
are  not  of  general  interest.  The  input  parameters  included  in  this 
section  are  for  the  most  part  shown  in  Fig.  1 and  are  listed  below. 


FREQL 

FREQH 


INCRM 


ISAME 


is  the  l_ower  frequency  limit 
is  the  h|gh  frequency  limit 

is  the  amount  by  which  the  frequency  is  incremented 
for  each  iteration  of  the  do-loop  which  chooses  the 
frequency 

is  set  equal  to  1 if  the  two  arrays  are  identical 
and  any  other  integer  if  they  are  not 


ER1,ER2,ER3  are  the  relative  dielectric  constants  of  the  three 
layers  of  dielectric. 


D1  ,D2,D3 


are  the  corresponding  thicknesses  of  these  layers. 
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All  frequencies  should  be  entered  in  gigahertz  and  thicknesses  in 
centimeters.  In  order  to  minimize  computing  time,  after  the  scan 
plane  is  chosen  (lines  59-61)  the  frequency  is  set  (line  67)  and  then 
calculations  are  done  for  each  angle  of  incidence  (chosen  in  lines 
84-93)  desired  before  the  frequency  is  incremented.  Some  of  the 
important  parameters  are  given  below. 

K is  an  identifier  representing  the  scan  plane 

(1  -+  E-plane  (<f>-plane)  and  2 -*  H-plane  (e-plane) 


LAMBDA 


B1  ,B2,B3 


is  the  frequency 

is  the  wavelength  in  free-space 

is  the  free-space  propagation  constant  3 


are  the  propagation  constants  61,62»3o 
Fig.  1 1 J 


shown  in 


PHII.THEI 


is  a dummy  variable  with  a value  equal  to  the  angle 
of  incidence  in  degrees 

are  the  angles  of  incidence  <t>i,0i  (in  free  space) 
shown  in  Fig.  1 (units  are  degrees) 


PHIIR, THEIR  are  the  angles  of  incidence  in  radians.  For  the  slot 
E-plane  (<t>-plane)  THEI  is  equal  to  zero,  and  for 
the  slot  H-plane  (e-plane)  PH  1 1 is  zero. 

The  second  section  starts  with  a list  of  input  parameters  (lines 
107-115)  which  describe  the  first  array  (the  one  between  dielectric 
layers  1 and  2).  The  parameters  which  specify  the  elenrnt  geometry 
are  shown  in  Fig.  17. 

L is  the  element  half-length 

W is  the  element  width 


T is  the  thickness  of  the  metal  sheet  from  which  the 

array  is  fabricated 

EFFRAD  is  the  effective  radius  of  the  element 

LTL  is  the  length  of  the  transmission  line  on  the  element 

(the  load) 
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ZTL  is  the  impedance  of  the  transmission  Vine 
DX,DZ  are  the  inter-element  spacings  shown  in  Fig.  1 


BDF  is  the  ratio  of  to  e. 

Again,  all  lengths  are  in  centimeters.  Lines  120  through  131  calculate 
the  total  admittance  of  the  first  array.  The  variables  introduced 
there  are: 

YS  the  self  admittance  of  the  array,  Ym(0)  [8] 

YI  the  admittance  of  the  images  involved  in  the  interface 

between  dielectric  layer  1 and  free  space  [8] 

YG  the  admittance  of  the  array  backed  by  dielectric  layer 

2 and  a ground  plane  at  distance  d2,  Y£  [8] 

YL  the  admittance  of  the  load  for  double-loaded  elements,  YL 

r 

YTl  the  total  admittance  of  array  1 > Y/(  + YL 

THEl.PHIl  are  the  angles  of  incidence  inside  dielectric  layer  1. 

It  should  be  noted  that  all  the  "admittances"  mentioned  above  have 
actually  been  multiplied  by  Z^/4  where  Z0  is  the  impedance  of  free 
space.  The  variables  which  begin  with  LE  and  DL  are  effective  element 
lengths  and  At's,  respectively,  where  L + DL  = LE  and  L is  the 
physical  length.  The  numbers  following  LE  and  DL  are  coded  so  that 
the  first  number  indicates  the  array  in  which  the  elenent  is  located, 
and  the  second  indicates  the  dielectric  layer  into  which  the  element 
is  assumed  to  be  radiating  for  the  particular  admittance  being 
calculated.  If  the  two  arrays  are  identical  lines  143-173  are  skipped. 
These  lines  calculate  the  total  admittance  of  the  second  array  and  are 
similar  to  lines  107-131.  Lines  178-184  set  certain  important  variables 
for  the  second  array  equal  to  those  for  the  first  array  in  the  case  of 
identical  arrays.  The  mutual  admittance  between  the  two  arrays  is 
computed  in  lines  188-190  where 

LEA2.DLA2  are  average  effective  lengths  and  At's 

YM  is  the  mutual  admittance  y|^  = Y^  . 

The  third  basic  section  (lines  191-288)  is  made  up  of  three 
subsections.  The  first  (lines  191-235)  finds  a total  admittance  for 
the  entire  sandwich  and  normalizes  it  to  yield  the  magnitude  of  the 
transmission  coefficient.  The  second  (lines  236-244)  determines  the 
phase  of  the  transmission  coefficient.  The  third  subsection  (lines 
245-286)  provides  for  printing  out  the  magnitude  and  phase  as  well 


as  plotting  the  magnitude.  In  the  first  subsection  lines  141-224 
determine  NORM  (line  224)  which  is  equavalent  to  the  normalizing 
constant 


s[*l  J h 

of  Eq.  (16).  Y (line  229)  is  equivalent  to  Y(d]  2 ,3 ;Yli  ,Y|_2)  in  Eq. 
(19).  Similarly,  TC(line  233)  is  the  numeric  transmission  coefficient, 
T (Eq.  (16)),  and  TDB(M)  (line  235)  is  the  transmission  coefficient 
in  dB  where  M identifies  the  angle  of  incidence.  The  second  subsection 
uses  subroutine  F to  find  F-j  and  F3  (lines  240-241).  The  phase  of  Fi 
and  F3,  FAZF1  and  FAZF3,  is  added  to  the  phase  of  Y(d]  o 3 i^LT »YL2) , 
FAZY,  to  obtain  the  total  phase  of  the  transmission  coefficient,  FAZT 
(line  244).  The  third  subsection  is  straightforward  formatted  output 
except  for  subroutine  GRAPH  which  is  not  included  here  because  it 
employes  several  plotting  subroutines  which  were  designed  specifically 
for  the  computer  system  for  which  this  program  was  written. 
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1 r ********************  wain  prograja  ************************ 

2 r * * 

A r * THIS  PROGRAM  CALCULAlrS  THE  TRANSMISSION  * 

H C * COEFFICIENT  a NO  PHASE  OELAV  POP  A BIPlANAR  ARRAY  * 

5 r * OP  Si  OTS  AS  S FUNCTION  OF  FREQUENCY.  ANGLE  OF  * 

6 r * INCIDENCE  AM"  SC  AIM  PLANE  A»E  ALSO  VARIABLE.  * 

7 r * * 

6 r **♦■♦♦#****’»************+****************************** **** 

9 r 

0 r EvTtRNALLY  COMPILED  SUBROUTINES  NEEDED 


INCLUDE  Sue.P9°9P;RCFLEC,29f °P I POISNB . 2Pr9P » 
2GRAPH«*?9E'>f,i 


PaRamETF RS  STORED  IN  COMMON  DATA  BLOCK 


ruMMON  PI.LAMpnA,aY.UZ«EFFROO,THEltPHlI.i  ,OXH,OZH 


REAL  LAM60A, L. I TL,  MAGIC,  TDB(9)  tLEll,LE12.LE23,LE??«LEA2 
RlAL  Yt(t>,2l'0),YH(*,200),  If'CPM 

roMPLEv  YS,  YI,vr.,YL,YM,YTl«  YT2,  y,F1,F3,  j.TC,R,ROOT,NORM,AA 
OOMPLEV  F l , F3 

DIMENSION  FAZYf'S)  ,FAZF1(5>  ,FAZF3(5),FAZT,a) 


DEFINE  ap*in  function 


akSini xy)=atan?(xx,sqpt(i .-xy*xx>  > 


OUTPUT  FILE  NAMES  AND  LOGICS 
NUMBERS 


rAl-L  ASSI  GN  ( 6Ho  ASHER  » 5H29B9M , ft ) 
faLL  AESlfcM*HTAH&A  ,5H29A9'VI.7> 


I'  PUT  DATA  AMD  CONSTANTS 


fkeql=* .0 

FREQHsl.2.0 

TNCRMsO.I 

NpOJNT  = ( FREOH-FRE  01.  )/INCPM 

PI  = 3.1'*159<'b3 

nn=PI/IB0. 

RD=180./HI 


SELECT  PLANE  OF  INCIDENCE 


i ^ ii  1 1 ■ .min min  i 


u 

I ’ 

I 

I 

l 

I 

I 


» 

I 


i ! 


it 


t 


56  r 

57  r 
*e  r 

59 

CU 

Cl  ° 
62  r 
62  r 

6*+  r 

65  P 

66 
67 
66 


ini 

102 

103 

104 

105 

106 
107 
106 

109 

110 


6 = 1 ....  Pul  plane 
y--d  ....  theta  plane 

K = 1 

f-0  TO  « 
k = P 

ITERATE  THROUGH  frequency  range 

nu  1 KKslf.NPOINT 
E = 1 

EkEOsF -!FwL  + KK*tNCR^ 

L AmB0A  = 50./FRfo 


69 

r 

70 

C 

COMPUTE  MrDlA  PROPAGATION 

71 

C 

72 

»s2«*PT /LAMBDA 

73 

»1=B*SQRT ( E R 1 ) 

74 

P2=B*SQRT  ( LR2 ) 

75 

R4=B*S0RT(EK3) 

76 

C 

77 

c 

SELECT  ANGLE  OF 

INCIDENCE 

76 

r 

M=1  ....  0 

t>EGRE£S 

79 

r 

=2  ....  20 

if 

80 

c 

=3  ....  60 

M 

81 

r 

=4  ....  70 

II 

82 

r 

=5  ....  60 

If 

83 

r 

84 

1 0 

GO  TO  ^2»3«4«5»6).w 

85 

? 

nuM= o.R 

86 

FO  TO  7 

87 

? 

nuM=30.0 

«6 

GO  TO  7 

89 

4 

PUM=60.0 

90 

GO  TO  7 

91 

* 

r>UM=70.0 

92 

GO  TO  7 

93 

6 

PUM=80.0 

94 

7 

TF(K,E^'.2)G0  TG  12 

95 

phII=OI'M 

96 

THPI=0.0 

97 

GO  TO  13 

98 

12 

PMlI=0. 0 

99 

THFI*0"M 

.00 

13 

thF1R=THE i*or 

PhTIR=phI1*OR 

OFTERNINF  AOMITT ANCF  OF  ARRAY  #1 

A p R A f PARAMETERS 

I =.375 
w=.ie 
Ts.0071 
fffrad=w/4» 


111  nx=1.3S5 

112  rz=1.3e>5 

113  I.  1 L s , 3 > 

11**  7Tl=24r. 

11b  RuF=1.26 

life  r 

117  r COMPUTE  El  EMENT  F.FFfcTIVE  LENGTH  Afn  ARP  AY 

116  r admittances 

119  c 

120  call  UFLL(L»W»T  » L AV80A  * OLll » FR1 ) 

121  IEll=L  + OLU 

122  TALL  POISOMFRI  ,1  ,EFFRAn, yS,LEll,DLll ) 

125  Call  HP  ISOM  fry »2»D1» YI ,LE11»PL11 ) 

12a  Call  0ELL(L,W,T, LAMBDA, DLl2,rR2) 

125  Lt  1 2=L  + DL  J.2 

126  CALL  POISoMFP2,5,02,  YG,LE12,DL12) 

127  Yl.  =0.b*J*ZTL*TAMB*BOF*LTL  )*RDF 
12b  C 

129  r SL'M  ADMITTANCES  FOR  TOTAL  ARRAY  ADMITTANCE 

130  r 

131  YTl=YS+YI+YG+YI. 

132  c 

133  r COMPUTE  ANGLES  OF  INCIDFNCE  IN  MEDIA  ttl 

13*+  C 

135  THFlrAoSlN (SIM (THEIR) /SORT ( ER1 ) ) 

136  PhI 1=AC  SIN ( SIM ( PHI J R ) /SQRT (ER1 ) ) 

137  IF<ISA"»E  .tC.DFO  TO  15 
136  r 

139  r DETERMINE  ADM  j TT  flMCr  OF  ARRAY  #2 

140  C 

141  r A ?R  AY  PAR  A.Mt  TrRS 

142  r 

143  l s • 375 

144  Ws.lS 

14b  T = . 0 0 7 1. 

14b  FFFRAD=W/ 

147  Px=1.3~5 

146  nz=l • 3fc5 

149  I TL=.3? 

150  7TL=24’, 

151  PUF=1.P6 
15?  r 

153  C COMPUTE  FI  EMEWT  EFFFCTIVE  LENGTH  Ann  ARRAY 

154  r AQMITTAMCES 

15b  C 

15b  TAI  L D«--LL  (LtvJ,T«LAMBDA,0L23,FR3) 

1*7  I E 23=L  + 0L23 

156  taLL  POISoN(FR3,1,FFFRAD, YS,I  F23,DL23) 

159  CALL  POISOM ER’ »2 » 03 » Yl »LE2S » Dl23 ) 

160  CALL  JELL(L«W,T, LAMBDA, DL22,FR2) 

161  » E?2=L+nL22 

162  CALL  H,**ISON(FP9,5,P2,YG«LE?2,0l22) 

163  YL=0.5«J*2Tl *TAM( B*BDF*LTL )*P0F 

164  r 

165  r SIM  ADMITTANCES  FOR  TOTAL  ARRAY  AO*"  T TT  AMCE 


I 
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YT2sYS+YI+Yb+Yl 


COMPUTE  AMGLES  OF  INCIOFNCE  IN  MEDIA  3 


ThE3=AR SIN ( SIN < THEIR > /SORT (ER3) ) 
PhI 3=APS1N< SIN (PhIIR) /SORT (ER3) ) 
rs(j  TO  1.4 


IF  ARRAY  #1  IS  ThE  SAME  AS  ARRAY  4?  THEM 
EQUATE  PARAMETERS 


LE22=LF12 

l E23=LF 11 

Pl_22=UL12 

Pl23s0i.11 

YTP.sYT] 

ThE3=ThE1 

PH I 3=PU  I 1 


DETERMINE  MUTUAL  ADMITTANCE  BtTFEEv  arrays 


l.fcA2=«LE12*LE2P)/2. 

PLA2=(i'La2  + L'L2p)/2. 

Call  PPISON(PRP,f ,P?. YM»LLA?,DeA2) 


COMPUTE  NnRMALl2ATI"\  FaCTOP 


7R (K.t0.2)bOTO  21 
ri;SFl  = COS(PhIl  > 
roSF3=CDS(phI3) 
r-o  TO  22 
PuSF1  = C0S<THE1  ) 

C0SF3=C0S(ThE3) 

AA=24U.*PI*CEXP(-J*fll*EFFRA0*C0SFi)/<B*«*0X*n?* 

1 RuRT ( C"SP1*C0SF3) ) 

PRs ( COS ( dl*DLl1 )*C0S<B1*L*SIM«  THEl) )-COS»Bl*Lfll) 
1-SIN(B1*UL11)*SIN<THE1>*SIN<P1*L*SIN<THM  ) ) )/ 
^(ERl**r!«25*SIN(Bl*LEH)*C0S(THLl)  } 
CC=(C0S(Ba*UL?^»*C0S(H3*L*SIM<  THF3) ) -COS ( B3*LE23 ) 

1 -SIN(B'*UL23»*SII‘M  THE3)*SIN(P3*L»SIN(THf  3)  ) )/ 
2<eR3**1.25*SIN(R3*LE23)*CuS<TH£3  ) ) 

P00T1  = 1.-<1./E"I  ) *SIiM  ( THF1R  ) **2 
F1=SORT(hooT1  ) 

P00T3=P00T1-(1 ./E  PI )*SIN(PHIIR )**? 

S3sSQRT(K00T3) 

Call  RF FL1(3,S1 tS3,R»EPlt2> 

PuOT=( 1.+R*CEXP(-J*2.*P1*01*P0SF1 ) > / ( 1 . -p*CEyP ( - j*2 .*B1 *01 
P*C0SF1 ) ) 

S*RTl=SOhT  (PE At.  (ROOT) ) 

PO0Tl=l.-(l./E°3)*SIN( THEIR)*** 

S1=SORT(RooT1) 

P00T3=P0UTl-(l./t»*)*SIN(PHTTR )**2 
SasSQR  r( ROOTS) 

PALL  Rf FLI (3tSi #S3*R»ER3«2) 


221 

222 

223 

22*4 

22b  r 
2?b  r 
227  c 
22b 
22* 

23t  r 

231  r 

232  r 

233 

234 
2 33 

236  c 

237  r 
230  r 
23* 

240 

241 

242 

243 

244 
24b 

246 

247 

240  23 
24*  ?4 

250 

251 

252  r 

253  r 

254  r 

255  r 

256  1b 

257 

250  19 

259 

260 

2bl  1 6 
262 

263 

264 
26b 

26b 

267 

26b  2b 
26*  1 

270 

271 

272  r 

273  r 

274  r 
27b  1 7 


P00Ts< 1 .♦K*CEXP<-J*2.*P3*U3*C0SF3> > / ( 1 . -R *CE *P ( - J*2 . *B3*03 
?*C0SF3> ) 

5uPT3=S0KT (REAL (ROOT) ) 

* qRH=1./( AA*RP*CC*SQRTl*SQrtT3> 

CPMPUlE  COMPLEX  ARRAY  ADMITTANCE 

rpHl22=B*D2*C0S(UUM*DR ) 

Y=(  vri*r  r<'-TM*YM)*CEXP(-j*CPHl22)/YM 

compute  transmission  lOfff.  ano  convert  to  ob 

TC  = 1./»M0RM*Y  ) 

F AG  TC=r  APS ( T C ) 

TuB  < M ) =2o  #*ALO<~10  ( *-:AbTC  ) 

COMPETE  PUASE  DELAY 

PaZY(M)=ATAM?(  AIM  A|?(Y)  « REAL  ( Y ) ) *RD 

tall  P (Ol.K.TprJ .THEIRiPHI1.PHIIP.ER1iH.o1.fi  ) 

TALL  FfD3.K,THF3.THCIRtPWl3.e>HiIR«ER3.B,R3«F3> 

FaZFI < >=ATAN?( AIMAG(1./F1 ) ,REAl (1./P1J )*RD 
FA7F.M  )=ATAN2« AI«AG(1./F3) »PFAL(1,/F3) ) * R 0 
Fa2T(D)=FA^Y(M1+FAZF1((“  )+FAZr3(M) 

TP ( K • E ' .2)60  T n 23 
Yt ( M.KP  +1 ) =TOP  ( M ) 

PO  TO  24 

Ym(M,K- +1 )rTOP("> 

r=M+i 

Tp (M.E  .6 )G0  TO  1 6 

n o To  'o 

P(’  I TE  TRA^SM  JSSION  COEFFICIENT  AND  pHASF 
rFL/Y  OUTPUT  FTLES 

TP  (K.E1  »l.Ai\C»KK.EO.O)  WRITE  (6.1*) 

IF ( K • £'  .l.Ar.O.XK.tl'.OjP.RlTP  (7,i9) 

PuFMAT(  ?X, ‘SLOT  E-PLAIME  (PHI  PLANE)') 

TP(K.E-.2.  AMH.KK.EO.O  )WlRTTt  (6«lfl> 

TP  (K.Er  .2.A  '’0.KK.f  0,0) WRITE  (7,16) 

P JKMAT(2X , 'SLOT  F-PLANE ( ThETA  PLANE)') 

!v'HITE(‘-  ,*0)FPFO,  (TO0«M)  ,*=1,5) 
rib  I TE  ( 7,20>FPE  ''t  (PAZ  Y ( r ) ,|W=1,5) 

LkITE 17,26) (FA7F1 (4) ,w=l,s) 

'KITE (7,26) (FA7F3(- ) »M=l,b) 

P hITE(7, 26 ) (FA7T( ,Msl»5) 

P'hl  T E ( 7 , 2. 6 ) 

F tjRMA  T ( 1A  x , 1 (FR.2,33)) 

FONTIN'  E 

TP ( K.Ef  .2)60  TO  1 7 

r-o  TO  ° 

Cl  OSE  OUTPUT  PILES 


r( GSE  f 
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*7o 
277 
27B 
279 
2«0 
2H1 
2 St 
286 
26H 
266 
28b 
267 
286 


r COSE  t 
PaUSE 
r 

r E ECUTf.  PLOTTING  SUdROUTInE 

r 

Tall  C-  ■''IPHlFPFOL  f et}fc  ilH,  iNChM,  1.25.1.C.0.U  »S.p.1*1«-35.  , YF  *1,-1.  J 
tall  graph i f rfql  * freqh* Imcrm  « i . 25,1.0* c. 4 *5.0*1 *1 *-3s. , yh. 2,-1 . ) 
fc'hITE(r',*5) 

25  FORMAT  (2X,  »P|_OT  AGAIN  U=NO  1=YES») 

PtAO  >&*-) IPLOT 
TP( IPLOT.K&.1>rP  TO  17 

"6  r aLL  E>7T 

rwO 


1 


1 

4 

\S 


1 


- % 
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JU  mu  mj 
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B.  Subroutine  POISSON 


Subroutine  POISSON  computes  the  mutual  impedance  sums  which  are 
needed  to  determine  the  admittances  (actually  impedances  since  they 
are  multiplied  by  ZW4)  discussed  in  the  previous  section.  These 
mutual  impedance  sums  have  been  simplified  by  the  use  of  Poisson's 
Sum  Formula  as  discussed  in  [8],  and  a detailed  explanation  of  this 
subroutine  is  given  there.  A discussion  of  the  calling  parameters  is 
given  below  followed  by  a Fortran  listing  of  the  subroutine. 


ER  is  the  relative  dielectric  constant  of  the  dielectric  layer 
into  which  the  reference  element  is  radiating  for  whichever  impedance 
is  to  be  calculated. 


IMP  is  an  identifier  to  distinguish  between  the  various  impedances 
which  the  subroutine  can  calculate.  IMP=1  calculates  the  impedance  of 
a slot  array  immersed  in  an  infinite  slab  of  dielectric  with  a constant 
of  ER  (only  radiation  to  one  side  is  considered).  IMP=2  computes  the 
impedance  of  the  image  slot  arrays  caused  by  the  interface  between  free 
space  and  the  dielectric  slab  of  thickness  D and  dielectric  constant 
ER.  IMP=5  calculates  the  impedance  of  an  array  radiating  (to  one  side) 
into  a dielectric  slab  of  constant  ER  and  backed  by  a ground  plane  at 
a distance  D.  IMP=6  calculates  the  mutual  impedance  between  two  slot 
arrays  separated  by  a dielectric  slab  of  thickness  D and  constant  ER. 
The  other  values  which  IMP  can  take  are  discussed  in  [8]. 


D is  the  thickness  (in  cm)  of  the  dielectric  slab  into  which  the 
reference  element  is  radiating  except  when  calculating  the  self 
impedance  of  the  array.  For  this  case  D should  be  set  equal  to  the 
effective  radius  of  the  reference  element. 


Y is  the  output  of  the  subroutine,  i.e.,  the  desired  impedance. 


LE  and  DL  are,  respectively,  the  effective  length  and  the  dif- 
ference between  the  effective  length  and  the  physical  length  for  the 
reference  element  (in  cm).  Of  course  this  effective  length  is 
calculated  with  the  appropriate  dielectric  constant  (whichever  one 
the  element  is  radiating  into). 


I a 


I Ml 


coRROUTit  t POISSON  (ER»  IMP*  U»Y,LE«0«-) 

***+*******♦****•*■*♦****  *****************  *********************** 

This  suhRhutjmp  computes  all  op  thf  mutual 

TisPEOfl  'Ct  SL'PS  NEEDED  TO  UR  TFRmH''L  THE  FMTIRF  TER”INaL 
TmPlOA’CL  OR  OTPOLE  OP  SLOT  ARRAYS  COVERFO  L I TH  Cl- 
rtFCTR'C  SLiPS.  IT  WILL  ALSO  COMPUTE  ThF  SUMS  NEEDED  TO 
PlMO  T ME  TERMINAL  JMPFDANCE  o F A SLOT  ARpAY  COVERFD  yilTh 
A DIELECTRIC  SI  AS  PACKED  BY  A GROUND  PLANE « THE 
njTUAL  If  PE.OANCE  BETWEEN  TWO  Slot  ARRAYS,  AMI  the 
TfiPEDA1'  Ct  OR  A DIPOLE  ARRAY  SITUATED  ABOVE  A 

semi-infinite  dielectric  ground. 

these  mutual  impfoance  sums  are  computed 

p Y THL  USE  Of  POISSAN'S  SO.-i  pOrMULA  as  PfSCRTPEO 

in  the  accompanying  text. 

* IMPEDANCE  AMO  SCAN  FLAME  IDEMTIFJFRS  * 

IMP=1....  FOR  SELF  IMPEDANCE 

=?....  FOP  SLOTS  INSIDE  UlElECTRlC 
=*....  FOR  oiPOl ES  OUTSTOE  olELFCTRlC 
FOR  niPOlES  INSIDE  nrELECTRlC 
sF....  FOP  A GPOUivjD  PLANE 
= <=....  FOP  MUTUAL  AD|ViITT ANCF  BETLEEiV 
SLOT  ARRAYS 

Krl....  FOR  SLUT  E-PLANF  SCAN 
=?....  FOR  Slot  H-PLANE  scan 

******************************************************************* 

***  INrLUOE  EXTERNAL  RFFLECTTO.m  COEFR  . S.lBR. 


INCLUO-.  REFLEC  * ?9£9R 

***  PA'  Af'ETfRS  STORED  IM  C0M«"0N  DATA  BLOCK 
roMMON  PI ,LAMPOA.nx«D?»EFFKAn« fHFTA»PHI»l 
***  DEFINE  PARoML-TfR  TYPF 
PEAL  LAMBDA, LF,L 

roMpLf.v  J,SUM'1,SU«?,SUN3, TERM, SUM ,R,Y 
***  DEFINE  R UNCTIONS 


CoT(XX)=tOS(XX)/SI»  (XX) 

CoTH(Xv  ) = l./TA*'H<  XX  ) 

SINH(XY )=(tXP( vX)-EXP(-XY) )/?. 

»**  INITIALIZATION  AMD  PFPlNTTION  OF  PARaMETcRS 
***  FO'  The  VAPJOJS  cases 


.-sCO.t1  .) 

Ps2 « *P T /LAMBDA 
FHS=ER 
c E = 2.* 


5b 

57 

hb 

59 
6C 
61 
62 
63 

6H 

63 
6 b 

67 

68 

69 

70  r 

71  r 

72  r 

73 

74  1 
7b 

76  U6 

77 

76  r 
79  r 

60  C 
81 
Bid 
83 
P4 

8b  45 
66 

87  C 

88  r 

89  r 

90 

91 

92 

93  r 
99  r 
9b  r 

96 

97 
9£ 

99 

100 

101  r 

102  r 

103  f 
1 n 4 
10D 
106 

107  ?3 
inb 

109 

110  r 


TFT I^P.FO.l )DF=0 
CEl=Ot 

TF<IMp.FU,3)EPS=l  . 

IF  ( IMP.F<-.5>DF=EFPBAD 

TF  I I*ip« pu# 5 ) nri  =0  . i 

PtR=l./FHS 

S*EK=S -KT (t PS) 

RSRER=1  ,/SF.F  R 
pf=b*s(-ek 

SIi\|PhI=SHv(PHI*PI/160.  ) 

8lMTHE=SIM  THFTfl*PI/lsO.  ) 

r-bv,i=(''.«o. ) 

8:0^2=  ) 

90*13=  C'  . *0.  ) 

***  SE'RCH  FCP  POcctHLF  MI’S  (Mir  0 ♦ +1 « -1 « ♦ ?. -2 » . . • ) 

nu  13  ' M=0«2P0 
a i=NM 

i tf:sti=o 

roMTIiM'-'E 

PoOTl=T  .-PE.P*T‘?II';THE+Nl*LAWBnA/0Z>**2 

***  TtST  FCR  FOR  51  Rrfli 

TF (ROUTl.&T.O.TGO  TO  45 
ttFSTj=ITF.ST1  + i 

TF< ITEFTl.tP.PTPO  TO  34 
PU  TO  47 
roPTlWI  F 
9l=S«PT (K00T1 ) 

***  COMPUTE  RFAL  PATTERN  FACTOR 

PHF=(C  SlBt*OL>*COS<BF*L*RSRFR«<SINTHF+Ni *LAmB0A/02> ) 
?-COS(H«  *Lt  ) - F1M  (PF*OL  )*  (PSRF'R*  (SINTME+rJl  *L AfbOA/OZ  ) )* 
.^«lN('j£*L*PSnEP*(SI*;,TME  + Nl*LA^POA/[)Z)  ) ) /9i 

***  SEARCH  FOP  POSSIBLE  M3’S  <iV3=  0 1 + 1 * -1  t *? . -2  ♦ . . • ) 

no  24  f!W=0 « 200 
r'3=Nfi 
TTFST3=0 
roNTiro'E 

Pt)OT*  = 'OuTl-PF°*(SI mPhT  + M3*L ANrinA/OX  )**"?. 

**+  TEST  FCjW  S7  PEAL 

TF (R00T3.PT. n. >G0  TO  23 
TlR'6=(,,.«0.  > 

PO  TO  32 
S3=S*RT(R00T3) 

TF ( IhP . EO , 6 ) GO  TO  65 
TtPi*i=PRE*PFE*rrxP<-j*PF*nEi*ss>/S3 


***  c<J"PUTt  rvnOMt  NTIaL  TERM  FuR  VARIOUS  CASfS 


Po  TO(  i2»bl,4At54,f>7.60),IMP 
TALL  HI  FLU  ?«S‘»  »S$.R»EP«?) 

Tf.  R’ls  rj  Ri-i*P  / < 1 . -f’  * CfXP  ( -J*RF*Ofc.*S3>  ) 

PCI  TO  *?. 

Pfil.L  P^FLO(3,si  «S3,R.FR> 

rtP^=T'-  pr*p 

r,o  TO  32 

PALL  RFF|_I<3«SJ  »S3«R»FH»1  ) 

T£RW=T-  pf'i*R 
r,u  TO  -V 

7eRr=T-Rri*<-J*r0T(PE*0*S»>  ) 

Po  TO  *2 

TpRMr-  j*PRt*PRP/<  S3*SlMRt  *n*S3)  ) 


**«•  T£.CT  POP  COA'VfRGEMCF. 


TF  (CABS  (UR1'1  ) .i  T.  ( .O01*CftbS(«fUrtl)  ) ) ITFST  .s=lTrcT3  + 1 


***  CU  'POTE  FI^ST  SUMMATION 


SUM1=S"W1+Tfc  HP 

IF  ( IT6.FT3.E  0.2160  TO  ?° 

TF  ( N3.  IPO  TO  VU 

P'  — NM 

p-0  TO 

roWTifjot 


***  Rt'RCH  FOP  POSSIBLF  PR'S  (,*4=  0 , + 1 t - 1 -2 »...  ) 


T TFSTHrO 

COMTIMOF 

t>uOT4= pFK*  (SlNPHl+N4*L  Al^pyA/nx  )**£-RCOTl 


***  TtST  FC.F  Sa  Pt'AL 


TF (R0CT4.&T.0. IGO  TO  RR 

TfcRMst'..,  0.  ) 

00  TO  7 3 

S4=s.;ki  ( pocT 4 ) 

TFdITP.Ft.fc  J&O  TO  f-4 
Tt  P^=o»PHr.*PRF*FXP  I -HF*DF*S4)  /S4 


***  CO- PlJTE  EXPONENTIAL  TtRF'  FuR  VARIOUS  CASrR 


PO  T0(  «3.Si!.49,5b«SS»Sl ) ,IPP 
r aLL  RfFLI(‘*«si  .S4«KtER«2) 

TtF^sT- RNwR/ < 1 ,-F  *EXP(-PF*OF*Sh)  J 

pu  TO  *3 

PALL  RFFitM“«Sl  «S4,f?,rPI 
Tfc  "iRsT1  R i*h 
00  TO  '3 


166  ¥.5 

167 
16b 

16V  *6 
17U 

171  61 

172  C 

173  r 
17H  r 

17b  33 

176  r 

177  r 
17«  C 
17V 
160 
lai 
162 

163 

164  AO 
lftb  47 
166 

1?  7 

16b  13 
169  r 
19«j  r 

191  r 

192  *4 

193 

194 

19b  35 

196 

197  r 
196  r 

199  r 

200 
201 
202 

203  36 

204  r 
20b  r 

206  r 

207 
20b 
20<? 

210 

211  r 

212  r 

213  r 

214 

21b 

216 

217  39 
216 

219  r 

220  f 


CALL  KEFL1(«.S1 «S4,K.eH*1 > 

Tt^MsT*-  R^'*h 
PO  TO  33 

TfcPi»l=T»-  pi“.*C0TF<9E*ti*S4) 

Go  TO  33 

7t«l4=J*2K(:.*HRF./(S4*SlMH(RE*0*S4)  ) 


♦**  TEST  FOR  CONVERGENCE 


IF  (CABM  TERN). I T.  ( . 00l*C  ADS  ( SUfij. ) ) > lTEST4sITFST4.fi 


***  COMPUTE  SFrOP'D  SUMP" A T 1 ON 


Sum2  = S''32*TERM 

TF< ITEST4.E0.?)G0  TO  47 

T F ( N4  « L E • U ) G 0 TO  3C 

r.4=-Niv 

GO  TO  *1 

COP'TINUF 

IF (Nl.LE.O)GC  TO  13 

P'1S-NH 

GO  TO  46 

continue 


*♦*  St  WCh  FOP  POSSIBLE  Pi2  * S (..2  = 0 , ♦!  . -1  , ♦ ? , -2  , . . , ) 


TTEST2=0 
no  25  '■•1=0,200 
P 2 = nN 
CONTINUE 

RO0T2  = ^Ef'*(SIP'T9F  ♦N2*LAMRuA/nz )**2-1  . 


***  TEST  FCP  S?  REAL 


Tk  <ROOT2,(,T.O.  »SO  TO  36 
SUM=(0.  •(•«) 

GO  TO  37 
S*=SQKT(KOOT2) 


***  COMPUTE  IMA61  P'ART  PATTERN  FACTOR 


P l N = ( Cf  S ( BE  *0L  )*COS  (BE*L*KSRfR*(SINTHE+p,?*LAmBOA/UZ)  ) 
p-costb1  *Lf  )-STM(F.E*DL)*<RSRER*lSINTME+N2*LAM«DA/0Z)  )* 
aGiP'(ME*L*KS»ER*(Slf  THF<fN?*LAMBi)A/pZ)  ) )/R? 

9UMS(0. ,0. ) 


**♦  SLUM  FOP  POSSIBLE  PJb*S  (w5=  l) « +1 1 -1  • +?  « -2 « • • • ) 


TTESTo=U 
f>U  38  -,4=0,200 
P>b  = Nl4 
CONTINUE 

P()'1T5=^ER*  (SIP  DHl  + Nb^LAMROA/nX  )**2*R0CT? 


***  TEST  FOP  S6  PEAL 
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■ .! 


221 

r 

222 

2?3 

224 

22b 

40 

2?b 

2?7 

22b 

r 

229 

r 

230 

r 

231 

63 

23? 

53 

233 

234 

233 

50 

23b 

237 

23b 

6fe 

239 

240 

241 

6y 

242 

243 

62 

244 

r 

24b 

r 

24b 

r 

24  7 

4 J 

24c 

249 

26i 

251 

292 

293 

254 

36 

2“-b 

37 

25b 

257 

r 

25fc 

c 

259 

c 

260 

261 

262 

263 

264 

265 

?5 

266 

43 

267 

260 

269 

r 

270 

r 

271 

r 

272 

273 

274 

IF(R00T5.bT .0. )G0  TO  40 
tfPM=(1' . tc  . ) 

GO  TO  *1 
^3=SORT  (RpOTb) 

IF  ( I^.EG.hJGr  TO  63 
I t_RM=0*Pli“i*PIM*t>  P ( -h£*nF*Si5l  /S5 

»**  CO  ’PLTt  E *POl\<E>T IflL  TE.RP  FOR  VARIOUS  CaSfS 

r-o  T0(4i«b3*5r,b6,59»62>»i*P 
CALL  Ri Fill1 ,SPtS^,R,eF.9) 

TtP'tsT^  RF*R/(l.-P«E*P(-nr*Or*Sa)  ) 

GO  TO  11 1 

call  Ri-:FL0ts«sr>«S5.R*E:R) 

TlRWsT'-  R i“i  * R 
Go  TO  **1 

tall  RFFLI  (c  «SP»S‘,.R«FR»1 ) 

Tfc  Ri4  = T *•  Rt“i*R 

r,o  To  u l 

T£RM=TFRi*i»COTH  ( BE  *0*S5 ) 
r-0  TO  4 1 

TtPM  = J*Pi'^PIf'/(SD*SlOI  (PL*0*S3)  ) 

**♦  TL°T  FOR  CONVFRGtNCE 

TF  (CAUSdKRi1  ) .1  T.  ( .001*CABS<Rlli1)  ) ) ITESTk=ITf  ST5+1 

TF  (CABR(TE;hr  ) .FT.  i .001*C«BS(SU‘I)  ) .AMD.ITrST5.GT.O)ITFST5sO 

SuMsSU-’+TERi 

TF ( ITESTb.E0.4)GO  TO  37‘ 

TF  ( Nb  • C ET  • 0 • ) 00  TO  35 

r-o  TO  *9 
roWTiwME 

TF  <CABC  <SUF) .1 T. < .f  01*CARS<  SHM3) ) ) ITF ST?  = I Tf ST2+1 

TF  (CAB  - iSor  ) .FT.  ( .C01*CAPS(SliMd)  ) . AFD . I Tf ST? .GT . 0 T I TF ST2  = 0 

*•*  CO  PLTt  THT«l  SUf'i^ATlOiM 


S.U'rA  = Si  *5+SU* 

TF  (ITLOT^.t  ii,4lR0  TO  4? 

TF  (M2. 1 E .ii  KO  *0  ?t 
R'ir =— N»'» 

r-oTO  jb 
r o 1 » T 1 1 J ' . F 
roMTin  T 

ruilST  = 1?l!.**JI/(SRFR*(SIM(HE*l  F)  )**2*P»H*r>X*rj7  ) 
TF  ( IMP.Fw.?  )CO'>STrf  Of  ST*?. 


***  CO  'PUT t AP«IT T aMCF/T"PEPANCE  AND  KEThRN 


>sCOMST*(SUA  1 ♦SUf'  2-SU«’3  I 

PrTuKo 

rivo 
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C.  Reflection  Coefficient  Subroutine  Package 

This  package  of  subroutines  computes  the  reflection  coefficients 
for  plane, inhomogeneous  waves  incident  on  a planar  boundary  between 
free  space  and  a semi-infinite  dielectric  slab  of  dielectric  constant 
ER.  Subroutine  REFLO  (lines  1-52)  computes  the  reflection  coefficients 
for  an  array  of  dipoles  on  the  free  space  side  of  the  boundary,  i.e., 
outside  the  dielectric.  Subroutine  REFLI  (lines  56-110)  calculates 
the  reflection  coefficients  for  an  array  of  dipoles  or  slots  on  the 
dielectric  side  of  the  boundary,  i.e.,  inside  the  dielectric.  A 
detailed  explanation  of  this  package  is  given  in  [8]  and  a Fortran 
listing  is  included  below. 
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cU°ROUTir-.r  REF|  0<  MO«SA»SR.R«FR) 


Tl<I  S SUBROUTINE  cOMPUTFS  The  REFLECTION 

rueFFxntNT  on  the  fref  spacf  side  of  a ni- 

FlFCTRTC  SLAB. 


*********************** 

ro^PLEV  o.R.B 

J z ( 0 • < ? •) 

RFPsl./FK 
SHER=S  'RT (t  P ) 
RbRERs1 ./SLRT(FR) 
SAS=SA*SA 
Sb$=SB*SB 
P=SRER*Sfc+W*0. 


**  SEI  FCT  RFFl  FCTlOfj  COfFF.  SURSCRTPTS 


***  COMPUTE  »3i 


assort (1.-ffr*( i .-S9S> ) 

Ps(l./(SAS*(1.-SHS) ) )*(SBS*M . -RAS > * C A-R ) / ( A+P  ) + 
PJSAS-S-'SIM  Sb-SREP*A)  / ( SP+SRER*A)  ) 

PtTUR'J 


***  CO'  PUTE  R4 1 


A A = 1 • "PER* ( 1 .+RPS ) 

IF(AA.I  T.O. ) GO  TO  41 
AsSQRT ( AA  ) 

P=- J*B 
FU  TO  42 
AsSURT ( -AA  ) 

Rs(l./(SAS*(1.+SBS ) ) >*<-RBR*fl.-SAS)*<A-R)/( a+B)+ 
>»bAS  + Sr'S)*(RSRFR*a-SREP*A  )/(RSRER*B+SKE'5*A)  ) 
0£TUR(m 


***  COMPUTE  RSI 


AA  = 1.-PER*( 1,+eBS  ) 

TF( AA.LT.n. JGO  TO  51 

assort ( aa  j 

&s- j*b 

r,u  TO  * 2 

assort  < -aa  > 

Rs(-1./(SAS*< J ,+SRSI ) )*(-SBS*(l.+SAS>*( A-8)/(A+6)+ 
? (-SaS+SPS)*<KSpEK*P-SRFR*A)/(RSRER*B  + SRFR*A>  ) 

return 

F(\jn 


SUBROUTINE  REFLl < NO»SA»SRiR«fR»OORS) 

ir********-*********************************-****** 

this  subroutine  computes  The  reflection 

COEFFICIENT  CO*  ELECTRIC  OR  MAGNETIC  DIPOLES  IN 
MfcOIA  ? RflOIATTNG  INTO  MEOIA  1 

► ******:►**»***#*************************»******* 

ruMPLEX  J.R.A 

integer  uors 

Sr<rR  = StiRT  ( E K ) 

sas=sa*sa 

CUS=SB*SB 

P=SRER*SB 

***  select  reflection  coefficient  subscripts 

GO  T0(3t3t3«4»*')tNC 
***  CO'PUTt  H3? 


7b  3 AA=l.-t'R*<l.-SP$) 

77  TF(AA.LT.O.)GO  TO  31 

7«  AsSQRTl  AA)-K1*0. 

79  go  TO  3? 

60  31  As-J*Sc;RT  (-AA) 

Ml  *£•  IF (OOKS.LO.l )GO  TO  33 

fi2  f>=(1./<SAs*<1.-SBS)  J )*(SBS*<1.-SAS>*(A-RI/(A+H)  + 

63  ^ISAS-SPS)*{SH-GRER*A)/«SP+SRER*A) ) 

64  PtTuRN 

65  33  p=a./<SAS*(l.-SBS)  ) ) * < SPS*  ( I . -SAS  ) * ( SRfr*  a-SB  > / ( SREP*rt+SB  ) 

Gfc  ?*tSAS-SBS)*<l?-A>/0+A>> 

P7  RETURN 

80  C 

69  C ***  COMPUTE  Pu? 


91  4 A = SUP  T ( £K*  ( 3 ,+SBS )-l« ) ♦ J*u • 

Si;  TP  I UORS,  to.  I )&n  TP  41 

93  «=(1./(SAS*(1.+S'S) ) )*C-SBS*(l.-SAS)*(A-B)/(A+e)+ 

94  >(SAS  + S'<S)*(SP-SREP*A>/(SP4SRFR»A  J ) 

93  Rp  TURN 

9fc  4 1 RsU./f  SPS*  ( 1 .4SbS)  ) »*(-SBS*(l  .-SAS)*(SpcR*a-SB)/<SRER*A4SB> 
97  ^♦(SAS+RPS)*(B-P)/»B*A)> 

9e  PpTURN 

99  C 

101)  C ***  COMPUTE  PR? 

101  c 

102  S fl=S0RT(EK*(1.4«‘BS)-l,  )4J»u. 

103  TP (UORS. to, 1 )Grt  TO  51 

104  Ps(-1./(SAS*(1  .♦SBS>  > )*(-sPS*(l.-ESAE)*»f-B)/(  A+B)4 

105  ?I-SaS4CPS)*(SP-SPER*A)/(SB*SPER*A) ) 

10b  PpTURii 

107  Ml  ‘ =<-l./<SfS*C1 .4S8f ) » )»<-SBS*(1.4SAS)M(SREP*A-S8)/<SPER*A4SB) 

108  -SAS  + Sf  S)  * IP-A  ) /(cJ+A  ) ) 

109  OtTuRM 


This  package  consists  of  three  subroutines  - DELL  (lines  10-38), 

■!  SICI  (lines  40-72)  and  F (lines  73-92).  Subroutine  DELL  calculates 

a Li  to  be  added  to  the  physical  length  of  an  element  to  compensate 
for  the  dielectric  layer  into  which  it  radiates.  In  other  words,  the 
physical  length  plus  the  Li  from  DELL  gives  an  effective  length  for  the 
element.  The  calling  parameters  of  this  subroutine  are 

m 

L,W,T  which  are  the  pnysical  length,  width,  and  thickness  of 
the  element  in  cm, 

LAMBDA  which  is  the  free-space  wavelength  in  cm, 


DL  which  is  the  Li  (in  cm)  sought,  i.e.,  the  output,  and 

ER  which  is  the  relative  dielectric  constant  of  the  media 

into  which  the  element  is  radiating. 

Subroutine  SICI  is  used  by  DELL  to  calculate  sine  and  cosine 
integrals . 

Subroutine  F calculates  the  F^  and  Fh  functions  (see  Eq.  (A-3) ) for 
slot  arrays  imbedded  in  a dielectric  slab.  The  calling  parameters  are 

D which  is  the  dielectric  thickness  in  cm 

K which  is  an  identifier  equal  to  1 for  E-plane  (<}>-plane) 

scan  and  equal  to  2 for  H-plane  (e-plane)  scan 

THED.THEI  are  theta  angles  of  incidence  in  the  dielectric  slab 
and  free  space,  respectively 

PH  I D -,PH  1 1 are  similarly  the  phi  angles  in  the  dielectric  and  free 
space 

ER  is  the  relative  dielectric  constant  of  the  slab 

B,BE  are  the  propagation  constants  in  free  space  and  in 

the  dielectric,  respectively 

FF  is  the  desired  F function. 

Note  that  the  angles  of  incidence  should  be  entered  in  radians. 


************  »»ism  LANEOUS  SUBROUTINE  PACKAGr 


THIS  PACKAGE  CONSISTS  OF  TwO  SUrPOUTtNES • * 

[JELL  AND  STCI  USED  BY  THF  MAIN  PROGRflM  * 

FOT  TmE  CAI  CULATION  OF  TuE  AOMITTAMCr  * 

OH  IMPEDANCE  pOR  arrays  nF  St,OTS  OR  oIPOlFS.  * 

* 

SUBROUTINE  DELI.  (L.w.T.LAmBPA.UL.ER) 

T*  IS  SURRn'JTl'E  CALCULATES  A LENGTH  INCREMENT 
TO  oK.  Appro  Tn  THE  PHYSTCAL  LENGTh  nF  THE  SLOT 
no  COMPENSATE  for  the  oIlLECTRIC  layer  INTO 
WHICH  IT  SATIATES)  TO  GTV£  A TOTAL  EFFECTIVE 
LENGTH. 

REAL  L.  LLAM.KHAT ,Kl*K2,LAMBPA,LEtLL 
COMPLEX  /m,2MM,2Ms:»200.ZA,ZAA.2l. J 
J= ( 0 . 0 , 1 . O ) 

P1  = 5.11*1S,3 
P=2.U*hI*6GKT<fR ) 

I LAhsL/LARBOA 
AM 

call  sici isibl .cibl*bd 

S 1 t'L=S  I 3L  *F  1 / ? . 0 

Call  siCi(SI?*T|  »CI2Bl»2.*BL  ) 

'•1?BL=SIcbl+PI/2. 

Call  Si  Cl (SI4PI  .CI4BL ,4.*bL > 

RI4BL=sI4HL+PT/2. 

kmAT=1  hi,  * ( -ALOGI  PI*W/l,  Aai«uA  ) +CIPL  + 0. <+i»2p  + 0 . 5*  ALOG  12.0  ) ) 
yHAT=fa!.*SI2HL+SO.*COS(2.*BL)*(2.*Sl?PL-Rl4PL)-3C.*(ALOG(BL/ 
24.  )•*■(!.  -.772-Cl4RL  + 2.*CI2BL-2.*ClBL  )*SIN<P.*Bl  ) 

PFL  = XHAT/(B*KHAT)  + ( l.ttt  -03)*W*t\HAT/(l  AI"!ppA*AlCG(3.*W/(2.*T 
2*  ) > 

pL  = OEL*LA!'"BPA 

PeTURivi 

E.gD 

SUBROUTINE  SICT(SI.CI*X) 


T i IS  SUhROUTIpE  IS  USEP  PY  DELL  TO  COMPUTE 
STNt  and  COSINE  INTEGRALS. 

7=ABSI v ) 

I F ( Z-4 . 0 ) 1 • 1 . 4 
v=t4.U-Z)*(4.0*Z) 

SI=-1.  -7U797E0 
IF ( Z ) <4  , ? » 3 
CI=-1 .0E3ft 

return 

ST  = /*I(  ( ( (1 .7*314  ir-9*Y  + l.  56  A9<JBF-7)*Y  + 1 .3741  68E-5  ) * Y 
r*h. 9394ft VE-4 ) *y*J .964B0?E-2)»Y+4.395S09c -l+ST/X ) 

CI=(  ( S . 7721 S6F.  -1  ♦ A l OG  ( Z ) )/Z-7*(  ( ( ( (1 . 3fcAqB5t-10*Y 
C*i.5ft4',Dht-A )*Y+I . 72S742r-b)*Y*l.l05999r_4 )*y*4.P90920E-3)*Y 
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V 

! 

— j 

5b 

01.31  5^0«f-l  ) >*Z 

57 

Pl.TUKi'J 

5B  <4 

fi=sin(z> 

59 

y=ros< 

• 

6D 

7=4.0/? 

61 

!=((((«((  ( 4.  04O069f-3*?-2. 2791431-  -2 > *7+5 . 515p70l-2 > *? 

62 

07.261.  4«'E-2>**+4.9e7716F-?)*?-3.3325l9F_3)*7-2.3146l7r-2)*2 

; 

65 

C-l .13+558E-5 >*7+6.75001 lr. 2) *z+2.5839B9r-10 

f. 

6*+ 

>/=(<((<((<< -5.1 06 69°£ -3*7  + 2. el9l79£-2 ) * 7-6 . 537283F-2  >*7 

65 

r+7.902T  34^-2) *7-9. 40041 6r-g) *Z-7.995556o _3) *7+2. 60 1293E-2 )*Z 

66 

0.1.764  OOF -4 ) 7-3. 12241HF-2 )*7-b.£9b441t-7)*7+2.500000£-l 

f 

r I=2*(  - I»V-Y»U> 

60 

51  = -?*  (SOl>+Y*V) 

69 

TF(X)5.b*6 

7 u 

Si=-3. 191593FP-SI 

71  6 

R£ TURN 

7<L 

Fn/D 

1 r 

this  surhouiinc  calculates  ff  and  fm  fop  slot  arrays 

2 r 

TM«tnOFn  IN  A niELFCTRlC  SLAP 

5 

COMPLEX  J.FF 

m 

roT(xx)=coS(yx)/ST''i(xx) 

1 

.|=(0.«1.) 

1 

b 

T6  CK.£C.2)GC  Tn  1 

n 

7 

PhI02  = -aTaA  (-SORT  ( fR  ) *C0S(P(4TI  )*roT(  f £*r  ,C0S  (f  HID  ) ) / 

I 

6 

proS(PHin) )+h*n*ros(PNii ) 

9 

FF  = -J«TEXP<  J*PuK  ?)  *SORT (£K  ) *C  v)S(PHTI  )*foS(F'*D*CP*-(PhTT  )-PHl.2)/ 

10 

2<rOS<PHTLl)*5INf  pf  *n*COS(PHlP)  n 

11 

Rp  TURN 

12  1 

Th£02=-AT  AN ( - SORT  ( f:R)*COS(  THrn)*roT(PE  *n*COS  ( THEN ) ) / 

13 

7rnS(THF  I > )+t'.*r>*COS(  Thtl) 

■ 

1*+ 

FF  = -J*f£XP(  J*TmF(»2)*SqPT  (t  H )*CuS(  THFD)*FnS(!<*n*CnS(  TH£T  )-lHE02>/ 

15 

2(rf)S(  TP-Fl  )*SIN»PF  *n*COS(  THr  r ) 1 ) 

: 

1& 

Pp  TURN 

1 7 

F ,\|13 
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